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ABSTRACT 




A VLF relative navigation system makes use of the fact that, at 
any given point on the earth, phase delay of a received VLF signal is 
highly stable and predictable. As the receiver is physically moved, 
phase delay changes linearly with distance from the transmitting sta- 
tion, so that by keeping track of the phase delay of the received 
signal from several VLF stations, one may keep an accurate plot of 
geographical position. 

Two problems experienced in measuring phase delay are the diurnal 
shift and long path interference. An investigation is made into a 
method of predicting the diurnal phase shift and the resultant phase 
due to simultaneous reception of long and short path signals. Also 
investigated is a receiving antenna having a cardioid shaped radia- 
tion pattern which could provide discrimination against long path 



signals . 
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1. Introduction 



Certain characteristics of the Very Low Frequency portion of the 
electromagnetic spectrum make feasible many important applications, some 
relatively new and others not so new. One striking application is the 
transmission of radio waves to submerged submarines which have not even 
an antenna above the surface. Another is the transmission of messages 
to mobile units operating halfway around the world. 

The application which inspired this paper, however, is that of a 
VLF relative navigation system. Such a system, using VLF transmitting 
stations already broadcasting for other purposes, has been tested and 
found feasible for both ships and aircraft operating at long distances 
from the transmitting stations. The system makes use of the fact that, 
at any given point on the earth, phase delay of a received VLF signal is 
highly stable and predictable. As the receiver is physically moved, 
phase delay changes linearly with distance from the transmitting sta- 
tion, so that by keeping track of the phase delay of the received signal 
from several VLF stations, one may keep an accurate plot of geographical 
location . • 

This paper attempts to resolve one particular problem associated 
with the relative navigation system application. This problem is the 
result of interference noted from the reception of a VLF signal which 
has travelled along a great circle path more than halfway around the 
circumference of the earth, i.e., interference of the so-called long 
path signal with the short path signal. Section 2 is a study of as- 
pects of propagation peculiar to the VLF spectrum, while Section 3 
discusses in detail the VLF relative navigation system. Sections 4 
and 5 provide information on VLF receivers and frequency standards 
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applicable to the navigation system. 

Discrimination against the long path signal is desirable to mini- 
mize navigational errors; this may be accomplished easily aboard ship 
by the use of an antenna system having a cardio id-type radiation pattern. 
VLF antennas are discussed in Section 6, with emphasis on the use of a 
loop in conjunction with a whip antenna to provide the desired cardioid 
pattern. 

Sections 7 and 8 provide a detailed discussion of the prediction 
of attenuation and phase delay of VLF waves. These predictions are 
necessary if one is to calculate in advance the effect of the long path 
signal on the observed resultant phase delay at any given geographical 
location. Section 9 provides an example of just such a calculation. 
Observed data on long path signals could not be obtained locally, since 
there were no VLF stations at a sufficient distance from Monterey, Cal- 
ifornia, to provide a measurable long path signal. It was thought 
that station GBR, transmitting on 16.0 Khz from Rugby, England, would 
provide this data; however, this station was off the air during the 
period available for experimental investigation (July -November 1966). 
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2. Characteristics of Very Low Frequency Radio Propagation 

The characteristics of VLF propagation which make these frequencies 
useful for long distance communication and highly accurate (to within a 
few parts in 10**) frequency comparison are a relatively low path atten- 
uation which is stable with time and a phase shift which is very constant, 
and predictable, over long periods of time. VLF also overcomes (to a 
large extent) the two major disadvantages of HF transmissions, unrelia- 
bility during ionospheric disturbances and fairly rapid and deep fading 
[28]. 

The major disadvantages of VLF transmissions are the very narrow 
(20-150 Hz) bandwidths available and the large transmitting installations 
that are necessary. The narrow bandwidths available are due to the small 
spectrum (10-30 kHz) and, mainly, to the high Q of VLF antenna systems. 
Since these antennas are very short with respect to a wavelength their 
input impedance consists of a large reactance and a very small radiation 
resistance. To obtain maximum efficiency, therefore, it is mandatory 
that all losses in the antenna, loading coils, and ground be kept to a 
minimum. The resulting low resistance coupled with high reactance gives 
these systems an inherent high Q and consequent narrow bandwidth. 

The high powered transmitters are required to obtain large areas 
of coverage and because of the high level of atmospheric noise present 
at VLF. This atmospheric noise is due mainly to lightning discharges 
which are centered in frequency at about 11 kHz. A typical discharge 
will radiate approximately 10^ joules and have a 3db bandwidth of about 
12 kHz [6j. 

Theory of VLF Propagation 

The two theories that have been predominantly used to explain the 
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characteristics of VLF propagation are the ray theory and the waveguide 
mode theory. It has been stated [ 28 ] that the mode and ray theories are 
related and, if all modes and all rays were taken into account with ap- 
propriate corrections for losses, both methods would yield the same re- 
sults. However, at long distances the mode theory is simpler to use 
because only one mode need be considered but many individual rays must 
be taken into consideration. At short distances many modes are present 
making the mode theory calculations complicated, but only a few rays need 
be considered since the higher order rays are incident on the ionosphere 
at much smaller angles where the reflection coefficient is small. In 
this paper the only distances which will be considered are those for 
which the mode theory calculations may be assumed valid with only one 
mode present (i.e., greater than 4000 km from the transmitter). 

From observed measurements of the electric field within the earth- 
ionosphere waveguide it has been found that the model shown in Figure 
2-1 can be used to approximate the actual physical case. The mode of 
principal interest at long distances is the TM^, shown in the left side 
of Figure 2-1, the TM 2 and higher order modes being small enough in 
amplitude beyond 4 Mm to be neglected. From Watt and Croghan [ 18 ] , for 
distances greater than 1 Mm the vertical electric field can be represen- 
ted as the sum of a very few waveguide modes: 
n = 1 to 3 



E z = Z E z ,n d>1 

n = 1 

Here E is the electric field strength for each particular mode (in 
z , n 

db relative to 1 v/m) 



E z [db, v/m] = 10 Log P r + - 10 Log f 
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Electrical f.ield structure for waveguide modes. 



(Wait CUD) 
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- 10 Log h i}t + 10 LogA njt 

- 10 Log (a sin d/a) - 10 Log h • 

L 9 r 

* 10 LosA n>r * 20 log G hjn 



- L, -txd/10 6 
d,n 

where 

P r = power radiated into the half space above the earth (watts). 

K*n = 10^-3 db for modes 1,2, a constant which relates field 

strengths at the surface to power radiated, 
f = frequency (Hz) 

h_- . = effective ionospheric height at the transmitter (m) . 

L , L 

h^ r = effective ionospheric height at the receiver (m). 

K, t = mode excitation modification factor, defined as the ratio 
of power launched into the concentric spherical shell guide relative to 
that for a flat guide with perfectly conducting earth, 
a = earth's radius (2^6.4 Mm) 
d = path length (m) 

(a sin d/a) accounts for energy spreading in the spherical 
guide structure. 



n,r 



= mode field modification factor for a curved earth at the 



ionospheric height at the receiver. 

20 Log G^ n = height gain factor (db) relative to the field at 
the surface. 



n = loss due to discontinuities along the path (db) . 

= effective attenuation rate in db/Mm (different for each 

mode) . 

These factors are discussed at some length by Watt and Croghan and will 
not be commented on further in this paper, with the exception of the 
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attenuation factor. 



For paths with uniform boundary surfaces, the attenuation factor is 
given by 

C< = M( (^ ) CXi + C< g 

where 

M( 0 ) = earth’s magnetic field factor 

= attenuation with magnetic field 
attenuation without magnetic field 

0 = direction of propagation with respect to the earth’s 
magnetic field 

= attenuation constant of the ionosphere 
CX g = attenuation constant of the ground, 
and CXg are related to the reflection coefficients of the earth and 
the ionosphere which are in turn related to the conductivities of these 
two media. 

The attenuation curves shown in Figures 2-2 and 2-3 are derived 
for M(0) = 1, i.e. independent of the earth’s magnetic field. Here 
H is the conductivity gradient of the ionosphere. From these curves the 
difference in attenuation of the first two modes is apparent. At 15 kHz 
and for a 4 Mm path the smallest difference in signal strength between 
the two modes is about 13 db, justifying the earlier assumption that 
only the first mode need be considered at these distances. Figures 
2-4, 2-5, and 2-6 are shown for comparison and while there are some dif- 
ferences in the attenuation curves for n = 1, the n = 2 curves show the 
same larger attenuation rates (with respect to the n = 1 curves). Fig- 
ure 2-7, while not a set of calculated curves, shows very vividly the 
effect of the earth’s magnetic field on the attenuation factor. 
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The phase velocity of a certain mode for a fixed set of conditions 
can be determined from a solution of the mode equation. This is done 
by Wait £ll3 and Wait & Spies resulting in the curves shown in 

Figures 2-8 through 2-11. The effect on phase velocity of several fac- 
tors can be seen here. Assuming a constant frequency, it can be seen 
that the height of the ionosphere has a large effect on the phase velocity, 
lowering it with increasing height. This is the cause of the familiar 
diurnal shift of phase delay of a received VLF signal. At night the D 
layer of the ionosphere (the lowest layer) dissipates and the effective 
height of the ionosphere is increased by about 20 km. The resultant de- 
crease in phase velocity causes an increase in the phase delay of the 
received signal as compared with a local standard. After sunrise the 
sun’s rays ionize the lower portions of the ionosphere, the D layer re- 
turns, and the process repeats itself. 

Comparing Figures 2-8 and 2-9 will show that an imperfectly re- 
flecting ionosphere (e.g., an exponential variation of ion density) will 
increase the phase velocity slightly. The effect of finite ground con- 
ductivity is seen in Figure 2-10, indicating that the phase velocity 
will be slightly lower over land (as compared to sea). From Figure 2-11 
it appears that the effect of the earth’s magnetic field on phase vel- 
ocity is small in comparison to the other factors. (The effect is much 
more pronounced for the n = 2 mode; however, that mode is being consid- 
ered negligible in amplitude here). Wait and Spies have many more curves 
of phase velocity, attenuation factors, excitation factors, etc., which 
would be of interest to anyone planning a more detailed study of VLF 
propagation. 
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Experimentally Obtained Data 

Figures 2-12 through 2-14 along with the experimental curves in 
Figures 2-2, 2-3 and 2-7, are examples of the vast amount of experi- 
mental data taken on the attenuation rates and phase velocities of VLF 
radio waves. The agreement of this data with calculated curves in cer- 
tain cases (e.g., Figures 2-13 and 2-14) is remarkable and provides a 
basis for using the calculated curves to predict phase velocities and 
attenuation rates. 

Figure 2-12 shows the many factors that must be taken into account 
in predicting attenuation rates: time of day, direction of propagation, 
frequency, and ground conductivity. In the 10-30 kHz frequency range, 
the range of attenuation rates is about 4 db/Mm (minimum to maximum) 
indicating that for long distances a considerable difference in signal 
strength could be expected depending on the path traveled. 

In this paper the experimental data will be used wherever possible 
and the calculated curves will be used as a check and to fill in the 
gaps not covered by experimental evidence (e.g., attenuation of the 
n = 2 mode). 

Summary 

The one characteristic of VLF radio transmissions of interest in 
this paper is the good phase stability of the signals over long periods 
of time and the predictability of this phase. Because of this phase 
stability frequency comparisons over long distances are possible with an 
accuracy on the order of one part in 10^. The amplitude and phase of 
these signals can be predicted fairly accurately at distances greater 
than one megameter from the transmitter by using the mode theory of 
propagation. The actual solution of the mode theory equations is a 
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rather complex problem; however, the solutions of Wait and others have 
been published and show a very good agreement with experimental evidence 
in some cases. 
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Figure 2-2 



Attenuation rates over predominantly sea water 
paths for n=l (experimental and theoretical). 

(Watt & Croghan [ 18] ) 




Figure 2-3 

Attenuation rates over predominantly sea water 
paths for.n=2 (experimental and theoretical). 

(Watt &. Croghan £183) 
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ATTENUATION, db/lOOOVm 




FREQUENCY, kc/s 

Figure 2-5 

Attenuation' rate of the second 
mode for a perfectly conducting 
earth and an imperfectly re- 
flecting iono sphere - 

(Wait [I *3) 
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ATTENUATION a, db/M^ 




Attenuation rate of* the first moo.e for an impel fee oly^ 
conducting earth and an imperfectly reflecting ionosphere - 

(Wait CUD) 




Figure 2-7 

Attenuation rates for 
predominantly sea water 
paths for moderate and 
low sunspot numbers (ex- 
perimental) - 

(Watt & Croghan C18J) 
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0.0 10 




FREQUENCY, kc/s 



Figure 2-9 

Phase velocity of the first mode 
for a perfectly conducting earth 
and an imperfectly reflecting 
ionosphere - 

(Wait CHI) 
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Phase velocity of the first mode for an imperfectly 
conducting earth an an* imperfectly reflecting iono- 
sphere - 



(Wait [ 111 ) 
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(v/c-l) x 100 




Figure 2-11 

Phase velocity as a function of frequency showing effect 
of reflecting height and magnetic field* 

'(Wait & Spies C 2 33) 
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Figure' 2-12 



.Attenuation rates deduced 
from various experimental 
data 

(Wait Cl3j) 




Figure 2rl3 




Figure 2-1 4 



Experimental phase velo- 
cities with pertinent 
theoretical data-. 

(Wait & Spies [23]) 
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Experimental phase velo- 
cities with pertinent 
theoretical data- 



(Wait & Spies G>;0) 
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3. A Proposed Very Low Frequency Relative Navigation System 

The proposed very low frequency (VLF) relative navigation system 
has been outlined by Stanbrough and Keilly [l6l, and presented in de- 
tail by Lake [253. The system has many advantages, providing (1) rela- 
tive simplicity; (2) a long range, accurate system with worldwide 
coverage; (3) the capability of using VLF phase-stabilized stations 
already on the air, with no requirement for additional stations; and 
(4) reliable 24-hour coverage. The only required equipment for ship- 
board use of this system are one or more phase tracking receivers and, 
preferably, a directional antenna system such as a loop used in con- 
junction with a sense antenna (whip). The latter is discussed elsewhere 
in this paper. One recorder channel to provide a continuous record of 
phase for each tracking receiver used would also be helpful. 

The extreme phase stability of VLF radio waves propagated over great 
distances, plus or minus one microsecond under normal conditions, is 
the foundation of this navigation system. An oscillator at the distant 
VLF transmitter provides the station with phase-stabilized emissions, 
while a stable oscillator aboard ship at the same frequency provides 
the local reference. This local oscillator is used in conjunction with 
a phase tracking receiver to provide an indication of relative phase of 
the VLF wave at the receiver. A high-quality quartz crystal oscillator 
would be suitable for use aboard ship and a typical unit in use at the 
U.S. Naval Postgraduate School provides a frequency drift (due to crystal 
aging) of about 2.5 parts in 10 ^ per day. This drift would be deter- 
mined for the particular oscillator in use aboard ship and its frequency 
would then be corrected on a regular basis, perhaps daily, to keep it 
within some specified frequency offset. At the time of taking a fix, 
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the phase variation due to the frequency offset of the oscillator would 
then be known and a compensation could be made. 

O 

Assuming a nominal value of 3 x 10 m/sec for the phase velocity 
of a VLF wave, the number of microseconds required for the wave to 
travel one nautical mile is found to be 6.18. In view of this, it is 
clear that for each increase in phase delay of 6.18^usec, the ship has 
moved one nautical mile farther from the transmitter. 

The VLF relative navigation system makes use of the foregoing; by 
keeping track of the phase of two or more stations as the ship moves 
away from an initial known point, an indication of the ship’s geograph- 
ical position is available at all times. By programming a computer to 
make the calculations and necessary corrections, one could even obtain 
a continuous track of the ship’s movement. It is, however, necessary 
to use VLF stations far enough away that mode interference is not ex- 
perienced, and only the first waveguide mode is involved (greater than 
4000 km) . 

One problem complicating the navigation system is that of the 
diurnal phase variation experienced with VLF waves . Fortunately, how- 
ever, this variation is predictable and the shipboard navigator can be 
provided with tables to correct for the phenomenon. Figure 3-1 illus- 
trates the diurnal variation, and is an actual plot of the phase of the 
13.6 kHz signal from the Omega Navigation System transmitter located at 
Haiku, Hawaii, referred to a rubidium vapor frequency standard at the 
U.S. Naval Postgraduate School. 

A first approximation to a solution of the problem of diurnal 
phase change may be made by assuming a trapezoidal pattern, or trapezium, 
as shown in Figure 3-2. A comparison with Figure 3-1 shows the similar- 
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Figure 3-1. Relative Phase of Haiku Omega Station, 
13.6 kHz, as Received at Monterey 
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Figure 3-2. Trapezium Approximation of Diurnal Shift, 
ity of the two patterns. Time intervals AB and EF have the same values 
of relative phase and represent the period of time during which the 
entire propagation path is in sunlight. Assuming the transmitter loca- 
tion west of the receiver, point B is the time of local (receiver) sunset, 
point C indicates sunset at the transmitter, point D local sunrise, and 
point E sunrise at the transmitter. It has been found that correction 
of sunrise and sunset times to those viewed from the effective height of 
the lower boundary of the D layer of the ionosphere (about 70 km during 
the day and 90 km at night) provides better correlation with the observed 
discontinuities in the phase plot. 

The diurnal shift (At in Figure 3-2) varies with distance from the 
transmitter, since it is caused by the fact that the phase velocity of 
a VLF wave has one value during the day and a lower value at night. As 
discussed elsewhere in this paper, these phase velocities may be pre- 
dicted for given paths, allowing the calculation of the diurnal shift 
at any given geographical position. 

Other, more precise, models of the phase plot may be determined 
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(as shown later in this paper). Whichever model is used, the procedure 
to be used in navigation is to predict the deviation of phase from the 
reference (all-daylight path) value, for the particular time of the 
desired fix, and subtract from indicated phase. As mentioned earlier, 
these deviations could be prepared and published in advance. 

An electronic navigation system such as that being considered is 
known as a circular grid system, as opposed to a hyperbolic system such 
as LORAN, which makes use of phase difference between a master and a 
slave station. In the VLF relative system, navigational charts may be 
overprinted with circles representing lines of constant time delay drawn 
about each VLF station as a center. These would be constructed using 
the 6.18 yusec/mile discussed earlier. (It should be noted that use of 
these charts introduces a small error since the phase velocity of a VLF 
wave is not exactly equal to the velocity of light; hence, actual phase 
delay differs slightly from 6 . 18 yUsec/mile. ) 

In order to use such an overprinted chart, one would first need to 
plot the ship's position from other information, which would then be 
used as a reference for future fixes obtained by the VLF relative sys- 
tem. The phase indicated by the tracking receivers for each VLF station 
being tracked would also be noted at the time the reference fix was ob- 
tained. At the time the VLF fix is desired, the phase readings on each 
VLF station would again be noted. For those stations whose propagation 
paths are all in daylight, the phase difference could then be transferred 
directly to the chart. This would be in the form of a time difference 
for each station which would indicate the distance travelled away from 
or toward the station. A short segment of a circle for each station 
would be interpolated on the chart between those segments already over- 



30 



printed. The intersection of these circular segments would then be the 
ship’s position. 

For those stations whose propagation paths are partially or totally 
in darkness, the correction for diurnal shift corresponding to the ship’s 
approximate position and time of day would have to be applied. There 
would be a correction for the time of the reference position and one for 
the time of the new position. The corrections would be subtracted in 
both cases, since phase delay is greater at night. After the diurnal 
correction has been applied, the phase information is used as described 
in the last paragraph. 

It is important to note that, since stations being used for navi- 
gation will be at least 4000 km distant, the circles of constant phase 
will appear as straight lines in small areas. In fact, if overprinted 
charts are not available, straight-line approximations of the circular 
segments may be used if the true bearings are known from the ship to 
the transmitting stations. These could be obtained from a previously- 
prepared table. The straight-line approximations will be adequate if 
the ship’s position is known within about 50 miles of its actual posi- 
tion [25]. 

A very complete example of a navigational fix obtained by this 
method, and using the trapezoidal approximation of diurnal phase shift, 
is given by Lake [25] . 

There are certain problems associated with the VLF relative navi- 
gation method, some of which are (1) cycle-slipping, (2) interference 
from the long path (around the world) signal, (3) inexactness of the 
trapezoidal model of diurnal phase change, (4) transmitter outages, 

(5) frequency offset of the transmitter master oscillator, and (6) 
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sudden ionospheric disturbances (SID’s). 

Cycle-slipping is a phenomenon in which the tracking receiver indi- 
cates and tracks one or more full cycles ahead or behind the true rela- 
tive phase of the VLF wave. This often occurs near sunrise or sunset 
when relative phase is changing rapidly. Frequently the diurnal phase 
shift will be more than one cycle, so that cycle-slipping is especially 
likely to occur. Cycle-slipping may also be caused by interference from 
waveguide modes higher than the first order, and by interference from the 
long path signal [27] . 

The long path signal may be thought of as a phasor, having some 
magnitude and phase relative to the short path signal. Together, the 
two signals form a resultant which constitutes the composite signal re- 
ceived by an omnidirectional antenna. Figure 3-3 shows the phasor addi- 
tion of the two signals, with = |eJ /0° representing the short path 
signal (taken with phase = 0° for convenience), E 2 = |E 2 | /8 ° the long 

path signal, and E^ = |E P [ /(f° the resultant. The notation E is used to 
indicate a phasor quantity. Since the long path signal may travel 




Figure 3-3. Phasor Diagram 

hundreds or thousands of miles farther than the short path signal, the 
angle 0 may take on any value between 0° and 360°. It is clear, then, 
that for any value of |E 2 | which is an appreciable fraction of |E]j , the 
resultant will vary considerably from the short path signal. Since such 
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a condition will occur under favorable circumstances, the need for a di- 
rectional antenna which will discriminate against the long path signal 
is evident. The long path signal is particularly likely to be appreciable 
near the antipode of the transmitter, as well as at other locations where 
the short path signal for any reason undergoes high attenuation concur- 
rent with low attenuation of the long path signal. 

The trapezoidal model of diurnal phase change, discussed earlier in 
this section, is a good approximation in many cases; however, there are 
also frequent instances in which it introduces excessive error. A more 
refined model of diurnal phase change, useful in all cases, is presented 
in a later section. 

Transmitter outages may occur on a scheduled basis for maintenance, 
in addition to those unscheduled outages due to equipment failure which 
invariably occur. During the interval of a transmitter outage, the 
indicated phase of a tracking receiver remains constant. Since the 
transmitter master oscillator normally continues running during the 
outage, the receiver will resume tracking on the correct cycle when VLF 
transmissions resume, if the receiver has remained geographically sta- 
tionary. If, however, the ship has moved a distance on the order of 
some multiple of a wavelength, the receiver will track the incorrect 
cycle (the same one as if it had remained stationary). The problem of 
determining the correct cycle is known as lane identification, and is 
minimized in other systems by utilizing lower frequency (longer wave- 
length) signals. Additional information, such as the ship’s dead 
reckoning track or positional information from other means, is necessary 
to identify the correct lane, or cycle. 

Frequency offset of the transmitter master oscillator of a VLF 
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station is occasionally made, for example to bring the transmitted fre- 
quency into coincidence with the United States Frequency Standard. These 
frequency corrections would occur infrequently and advance notice thereof 
could be given to users of a VLF relative navigation system. With such 
advance notice, the navigator would have no difficulty applying a cor- 
rection for the offset. 

Sudden ionospheric disturbances are caused by three phenomena: solar 

flares, magnetic storms, and high altitude nuclear bursts. While these 
disturbances occur infrequently, they cause phase disturbances on VLF 
circuits over wide areas of the earth, causing errors in VLF navigation 
in the case of severe disturbances. Chilton, Steele, and Crombie [19] 
have pointed out that solar flares cause a depression in the effective 
height of the ionosphere, which in turn causes increased phase velocity, 
resulting in a lower value of relative phase at the receiving point. A 
typical flare has an onset time of one to 30 minutes and requires from 
one to two or three hours to return to normal. A typical phase change 
caused by the flare might be 10 microseconds, but varies with path 
length, frequency, and amount of ionosphere height change. 
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4. VLF Receivers. 



The receiver used to obtain experimental data on phase and amplitude 
of various VLF stations for this paper was a Tracor Model 599G VLF phase 
tracking receiver. A simplified block diagram showing the essential ele- 
ments of such a tracking receiver is seen in Figure 4-1. The output of 
the block marked VLF Receiver is at an intermediate frequency of 1 kHz 
with a certain phase 0 vlf’ an< ^ one ^ n P ut °f the ph ase comparator. 

The other input of the phase comparator, a 1 kHz signal of phase 0 A , 
is derived from the external frequency standard. The output of the phase 
comparator provides the electronic phase shifter with an error signal 
A 0 , which tells the phase shifter how much to shift the phase of the 
external standard to maintain phase tracking. The output of the phase 
shifter then goes to the frequency synthesizer which provides the ref- 
erence signal 0 A . The condition for phase tracking is A 0 = - 

"a = 

The phase tracking servo and long integration time of the receiver 
provide the capability of tracking with an input s ignal-to-noise ratio 
of -50 db. The rated sensitivity is 0.01 microvolts of signal at the 
antenna terminals to provide phase tracking. Using a loop antenna made 
by Tracor especially for this receiver (Model 599-600), this corresponds 
to a signal strength of 0.3 microvo Its /meter . With a 150 sec servo time 
constant (front panel switch provides time constants of 5, 15, 50, and 
150 sec) the equivalent noise bandwidth is 0.002 Hz. A front-panel 
switch on the receiver allows selection of either a broadband RF filter 
or any one of five narrowband ones. The narrowband filters provide an 
RF bandwidth of 500 Hz, while the IF bandwidth is 50 Hz. The receiver 
is tunable from 8.0 to 31.9 kHz in 0.1 kHz steps. 
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Figure 4-1. VLF Tracking Receiver 



The receiver has available several outputs that proved useful in 
experimental work for this paper: (a) 100 ^sec and 10 ^isec full scale 
recorder outputs, (b) a front panel bi-directional digital counter for 
cumulative phase change, (c) a 100 kHz coherent output, and (d) an AGC- 
derived relative carrier level recorder output. The 100 yusec and 10 
jjtsec full scale recorder outputs use the receiver’s internal phase com- 
parator to give a dc voltage level proportional to the phase difference 
between the external frequency standard and the received VLF signal. 
(Note: the phase comparator of the receiver does not actually operate 
at 1 kHz as indicated in the discussion of Figure 4-1. The receiver 
contains two linear phase comparators; one operates at 100 kHz to pro- 
vide the 10 ^tsec full scale recorder output, while the other operates at 
10 kHz for the 100 yUsec full scale recorder output.) These recorder 
outputs were convenient and proved quite satisfactory for continuous 
monitoring of the phase of any phase-stable VLF station. Moseley Model 
51^)0B strip chart recorders were used for monitoring these and the 
other receiver outputs, as desired. 

The electronic phase shifter of the receiver actually operates at 
10 MHz. One cycle of phase shift then corresponds to 0 . 1 yu.sec drift 
which feeds the front panel bi-directional phase counter. This digital 
counter provides a cumulative indication, so that by taking readings at 
two different times, one may determine by subtraction the exact phase 
drift that has occurred during the elapsed time. 

The 100 kHz coherent output of the receiver provides a 100 kHz 
signal from the frequency synthesizer section shown in Figure 4-1. This 
signal thus has a phase corresponding to that of the received VLF signal. 
This output proved useful as an input to an external phase comparator, 
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allowing a check on the phase of any local frequency standard (the other 
input to the comparator) with respect to the received VLF signal. Both 
the rubidium vapor frequency standard and the quartz crystal oscillator 
in use at the Naval Postgraduate School have 100 kHz outputs suitable 
for this (see Frequency Standards section). 

The relative carrier level output of the tracking receiver provides 
a dc level proportional to the AGC voltage developed in the receiver and 
is thus a relative indication of received signal strength of the VLF 
station. The AGC has a range of approximately 40 db, and peaks at about 
40 db above the minimum signal necessary for tracking. 

One problem occasionally encountered in experimental work was the 
very long integration time of the receiver. For example, a cardioid 
antenna system must be adjusted for proper phase and amplitude of the 
whip antenna when changing to a new frequency. To get a good indication 
of phase and amplitude, it was usually necessary to wait five to 10 min- 
utes after each adjustment of the whip. Since several adjustments are 
normally required, the process can be quite time-consuming. 

Adjustments are particularly time-consuming if one is receiving 
one of the Omega stations. Each of these stations transmits only about 
one second out of 10, so the long integration time problem is compounded. 
At the time of the experimental work for this paper, there were four 
Omega stations in operation (there are planned ultimately to be eight). 
Each station operates on a time-shared basis on two frequencies, 10.2 
and 13.6 kHz. They transmit an unmodulated waveform in a prescribed 
sequence for a prescribed length of time (e.g., 0.9 sec, 1.1 sec, etc.) 
to facilitate identification. Still in the planning stage for the sys- 
tem are the addition of 11.33 kHz by all stations on the time-shared 
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basis, as well as transmissions of 12.00 and 12.75 kHz (not on a time- 
shared basis) by some of the stations. In order to obtain phase and 
amplitude information on an individual station, it is necessary to 
clamp the AGC and blank the receiver during the time the desired station 
is not transmitting. In the laboratory at the Naval Postgraduate School, 
this was accomplished by a Tracor Omega Gating Unit. This unit may be 
used with four tracking receivers simultaneously, and a modification for 
the receiver is provided which allows switching between normal receiver 
operation and Omega operation. The unit proved highly satisfactory and 
much valuable data was obtained using Omega stations. 

It was useful in the experimental work to have an auxiliary VLF 
receiver available which could conveniently be tuned from station to 
station for quick audio identification. The receiver used was a model 
AN/SRR-11, and a similar type would no doubt be useful for auxiliary 
purposes in the shipboard VLF relative navigation system. 
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5. Frequency Standards 

A good frequency standard must be used at both the phase-stabilized 
VLF station and the ship using the VLF relative navigation system; hence 
it is a vital part of the system. Its use in the system aboard ship is 
shown in Figure 4-1. 

There are two basic types of precise frequency standards: the quartz 
crystal oscillator and the atomic frequency standard. The former is 
characterized by a predictable frequency drift as the crystal ages, 
while the latter has either an extremely small drift or none at all. 

The past twenty years have seen great achievements in the field 
of precise frequency control; today, time difference or frequency may 
be measured more precisely than any other physical quantity. Using an 
atomic hydrogen maser, a frequency stability (RMS deviation from the 
mean) of 2 parts in lO 1 ^, over a one-day averaging period, may now be 
obtained [30] . One problem which has arisen in the growth of technol- 
ogy of frequency standards is a lack of standardization of the defini- 
tions of certain terms. Two of these are short-term stability and 
long-term stability. The former refers to average frequency deviation 
measured over a time interval sufficiently short that long-term effects 
may be neglected. Since the averaging period is not standardized , it 
must be specified; it varies from a fraction of a second to perhaps one 
minute. Long-term stability refers to slow changes in average frequency, 
and is usually specified for crystal oscillators by stating fractional 
frequency deviation per day, e.g. , 5 parts in lO^-l per day. Since 
atomic frequency standards have little or no systematic drift, specifi- 
cations for them will often be given in terms of standard deviation or 
RMS deviation from the mean over a specified time interval such as one 
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day. 

Quartz Crystal Oscillators 

There are many factors which affect the frequency of quartz crystal 
oscillators, which must be taken into account or compensated for. In 
particular, frequency stability is a function of: (a) temperature; (b) 
time (aging); (c) stress, vibration, and acceleration; (d) drive level, 
and (e) nuclear effects. These effects are discussed by Gerber and 
Sykes [29]. 

Ambient temperature changes account for the principal frequency 
change effects in most crystal units. Current practice in precision 
oscillators is the use of a double proportional oven to control the 
temperature of the crystal unit. This provides a method of removing 
heat from around the crystal at the same rate it is generated# A sim- 
pler, but less desirable, method uses a thermostatically-controlled 
oven. In this method, temperature varies between a lower and an upper 
limit, whereas it is maintained constant when using proportional con- 
trol . 

The effect of time (aging) in a quartz crystal oscillator is easily 
observed by the variation in frequency deviation from day to day- Accord- 
ing to Berber and Sykes [29] , aging of thickness shear crystals is caused 
by four factors: (a) temperature gradient effects lasting several min- 
utes to several hours after a thermal disturbance, (b) stress relief 
effects, (c ) change of mass effects, and (d) structural changes due to 
imperfections in the crystal lattice. In all crystal oscillators there 
seems to exist an initial stabilization period, normally lasting a few 
weeks, in which the oscillator has not yet attained its ultimate fre- 
quency stability. 
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The effects on frequency stability of stress, vibration, accelera- 
tion, etc., may be reduced by careful design of mounting structures. 
Austin [24] has made an experimental study of the effects of vibration, 
shock, and swinging motion on frequency stability of several commercial 
crystal oscillator units. He concluded that stability of these units 
was degraded negligibly under conditions similar to those which might be 
encountered in a shipboard environment. 

Another factor affecting crystal stability is drive level, or ampli- 
tude of vibration. A typical 2.5-MHz crystal unit might exhibit a change 
in frequency deviation of from 1 part in 10^ to about 8 parts in 10^ with 
a change in drive level from 3 to 1000 microwatts [29] . For this reason, 
automatic gain control is used in precision crystal oscillators. 

Nuclear radiation also has an effect on the frequency of crystal 
oscillators; the extent of this effect is discussed in the literature 
[13] , [15] , and [2 6] . 

Figure 5-1, from Gerber and Sykes [29] , shows a schematic diagram 
of a modern precision crystal oscillator. The double oven is shown, and 
proportional control would be used. To maintain a constant drive level, 
a large amount of negative feedback would be used in the amplifier, in 
addition to the AGC. 

The Sulzer Model D2.5 is a crystal oscillator of the type shown in 
Figure 5-1 and was in use at the Naval Postgraduate School for the ex- 
perimental work of this paper. Its drift rate was determined to be 
about 2.5 parts in lO^ per day. Other units of this general type often 
have drift rates of about 1 part in 10^ per day, while similar units 
used to stabilize transmissions of some U.S. Navy VLF transmitters have 
average daily drift rates as low as a few parts in 10^ [29] . This 







Figure 5-1. Precision Crystal Oscillator 
general type of oscillator should be well suited for shipboard use 
in a VLF relative navigation system. 

Atomic Frequency Standards 

It is interesting at this point to note the difference between 
a primary and a secondary frequency standard. The former is a self- 
calibrating device which is intrinsically accurate. The latter, 
however, requires calibration against a known precise frequency 
before it can be considered a standard. A quartz crystal oscilla- 
tor may thus be classified only as a secondary frequency standard. 
There are several types of atomic frequency standards, including the 
hydrogen maser, the cesium beam standard, and the rubidium vapor 
standard. The latter, considered a secondary frequency standard 
since its frequency initially must be compared with some other 
standard, was used at the Naval Postgraduate School as an adjunct 
to the preparation of this paper. The hydrogen maser and the cesium 
beam standard are both primary frequency standards. All three of 
these atomic standards are now commercially available. 
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Atomic standards in general derive their precise frequencies from 
frequencies corresponding to the energy involved in the transition of 
electrons between energy states in particular types of atoms, such as 
hydrogen, cesium 133, and rubidium 87. The cesium and rubidium standards 
operate on the same basic principles; both have found wide use in VLF 
applications as well as other uses. Space does not permit detailed des- 
criptions, but since the rubidium standard was used in some of the ex- 
perimental work for this paper, its basic operation will be described. 

The heart of the rubidium standard is a high-quality quartz oscilla- 
tor, whose frequency is stabilized by a passive resonance cell filled 
with rubidium and an inert buffer gas. The operation of the resonance 
cell is based on an optical pumping principle. The rubidium vapor in 
the cell is illuminated by filtered light from a rubidium lamp. A 
microwave signal is applied to a microwave cavity surrounding the reson- 
ance cell. This signal is obtained by multiplying the frequency of the 
quartz oscillator. When the frequency corresponds to atomic resonance, 
energy transitions are induced in the rubidium atoms and a photo detec- 
tor observes these transitions through changes in the amount of light 
absorbed as a function of applied microwave frequency. A servo loop 
connects the detector output and oscillator, locking the oscillator to 
a frequency corresponding to the center of the resonance line. Since 
the resonance frequency is influenced to some extent by the pressure of 
the buffer gas in the resonance cell, this pressure must be properly 
adjusted when the standard is manufactured. The cell in then sealed 
and the unit is ready to serve as a frequency standard. 

Rubidium vapor frequency standards have been designed for U.S. Army 
tactical applications and for guided missiles and aircraft. They should 



also be readily adaptable to a shipboard environment. State-of-the-art 
characteristics for these standards, according to McCoubrey [30] , in- 
clude an RMS deviation from the mean of 5 parts in 10 12 for a one-day 
averaging period and a systematic drift of less than 3 X 10 per month. 

Cesium beam frequency standards are presently being used as the 
United States Frequency Standard (maintained by the National Bureau of 
Standards, Boulder, Colorado) and to stabilize the VLF transmitters of 
the Omega navigation system stations. While the short-term stability 
of cesium beam standards is not as good as that of rubidium standards, 
their principal advantages are excellent long-term stability and high 
intrinsic reproducibility. McCoubrey [30>] gives the state-of-the-art 

*1 o 

long-term stability as 2 X 10 ° per day (RMS deviation from the mean). 

Stations of the Omega navigation system effect an intercomparison of 

their cesium beam standards periodically (currently every four hours), 

- 1 1 

and are thus able to control their stability to better than 1 X 10 
per day [32] . 

The Varian Model V-4700 is the rubidium vapor standard that was 
used at the Naval Postgraduate School for the authors T experimental 
work. Over a period of about four months, this unit exhibited an RMS 
deviation from the mean of 1 . 36 X 10” 11 , using an averaging period of 
one day. 
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6. Antennas 



As mentioned in Section 3, the long path signal at times can intro- 
duce serious error in the VLF relative navigation system. For this 
reason, it is desirable to use some means of discriminating against this 
component of a received VLF signal. The cardioid antenna system, which 
has been widely used in radio direction finding and radio compass appli- 
cations, provides a ready solution to the problem. The basic theory of 
this system has been described by Terman L3l and others, and will be re- 
iterated and elaborated upon in this section. The loop antenna, basic 
to a cardioid set-up, will first be discussed, then the cardioid system 
itself, and finally the results of measurements on antennas at the Naval 
Postgraduate School will be presented. 

The Loop Antenna 

The loop antenna, in various forms, has been used for many years. 
One of its earliest applications was reception of stations in the broad- 
cast band when this service was in its infancy. Loops used for VLF 
work, as well as most other applications, have a diameter which is small 
compared to a wavelength. A loop may be circular, rectangular, tri- 
angular, diamond-shaped, or some other similar configuration. Kraus [2] 
states that far fields of loop antennas of different shapes are the same 
if the areas enclosed by the loops are the same. 

A vertical loop antenna produces a polar radiation pattern in the 
familiar figure-eight shape, as shown below in Figure 6-1. The xy-plane 
is parallel to the surface of the earth. Maximum response, of course, 
occurs when the plane of the loop is oriented parallel to the direction 
of arrival of the electromagnetic wave. 
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Figure 6-1. Loop Radiation Pattern 

Appendix III describes the theory of operation of the loop antenna. 
Equations are derived there for the relative response and phase shift as 
the antenna is rotated from 0° clockwise to 360°. Figures 6-2 and 6-3 
show plots of these equations for 0° to 180°. The frequency is taken 
to be 20 kHz, and the loop is assumed square, with the length of a side 
equal to 1.635 m. (This length was chosen because it was the measured 
length of a loop available for experimental work.) The loop is assumed 
parallel to the y-axis of Figure 6-1 and direction of arrival is from 
right to left along the x-axis as shown. The angle of rotation Oi is 
taken as the angle between the normal to the loop and the positive x-axis 
(this is shown in Figure III-2) , so that (X = 0° and 180° corresponds to 
zero response, while & = 90° and 270° are the points of maximum response. 

The graph of Figure 6-2 was arbitrarily terminated at -35 db for 
convenience, since output theoretically goes to - oo db at OL = 0° and 
180°. For rotation from 180° to 360°, the amplitude pattern is exactly 
the same as for the first 180° of rotation (this is obvious from in- 
spection of Figure 6-1). 

For the phase pattern, the second 180° was not plotted in order to 
show greater detail. At 0C = 180° and 0C = 360°, there is a discon- 
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Figure. 6-2. Loop Relative Output 
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FigUre.6-3. Loop Relative Phase Shift 
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tinuity of 180°. The pattern for Of = 180° to 01 = 360° looks the same 
as Figure 6-3, except it appears inverted, i.e., at 01 = 180°, — ► 

-180°; at 0L = -270°, Ip = -(180° + 3.4 x 10" 4 deg); and at OL = 360°, 

( Q — -180°. (The symbol — is used to imply the existence of the dis- 
continuities at OL = 180° and G£ = 360°.) The 180° phase reversal on 
the back side of the loop antenna is significant, for one must be care- 
ful not to rotate the loop antenna to the wrong side of the null if 
phase measurements are important. 

The Cardioid Antenna System 

A cardioid antenna system consists of a loop antenna and a whip 
(or ’’sense”) antenna whose outputs are appropriately combined to pro- 
duce a cardio id-shaped radiation pattern. Figure 6-4 is a simplified 
schematic of such a system. Both the loop diameter and the height of 
the whip are small with respect to a wavelength. Some provisions for 
tuning, which are not shown, would normally be installed at the base of 
the loop and at the base of the whip. There would also be provisions 
for varying the amplitude of the loop signal or the whip signal, or both, 
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as well as provisions for varying the whip phase a small amount above and 
below the 90° phase shift shown. These provisions are necessary to allow 
adjustment to obtain the cardioid pattern. The adjustment consists of 
making the whip antenna component of the VLF signal equal in magnitude 
and phase to the loop contribution, for the case in which the plane of 
the loop is oriented toward the distant transmitting station. 

Appendix IV contains the equations necessary to plot the phase and 
magnitude of the cardioid antenna system as the loop is rotated about a 
vertical axis. The necessary phasor diagrams are shown, and Figure IV-3 
shows the phasor diagram for proper cardioid adjustment, with the loop 
oriented toward the station. Proper operation of the cardioid system is 
independent of the location of the whip with respect to the loop, so long 
as it is not an appreciable fraction of a wavelength away. 

The radiation pattern of a properly-adjusted cardioid antenna sys- 
tem is sketched in Figure 6-5. Maximum response occurs with the plane 
of the loop antenna oriented parallel to the x-axis. 

DIRE CTION OF 
SIGNAL ARRWAL 

~X 

Figure 6-5. Cardioid Radiation Pattern 
Figures 6-6 and 6-7 are plots of the radiation pattern in Car- 
tesian coordinates, using equations derived in Appendix IV. The plot 
is for a frequency of 20 kHz and a square loop 1.635 m on a side. Ro- 
tation angle (X is the angle between the normal to the loop and the 
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Figure 6-6» Cardioid Relative Output (0°\0L<(l80 o ) 



52 




53 



(db) 



positive x-axis. The maximum response is 6 db, which merely indicates 
that the voltage delivered to the receiver is twice that from the loop 
or whip alone. The minimum response, with the loop pointed away from 
the direction of signal arrival, is approximately -63 db. In the forward 
direction, 6 db beamwidth is 180°, providing the desirable characteristic 
of a broad front lobe. 

It should be noted at this point that -63 db is not the minimum 
response obtainable on the "back” side of a cardioid antenna system. 

The ideal minimum response approaches -oOdb when the system is adjusted 
for this. The adjustment required for this entails adjusting the whip 
phase for minimum system response, which corresponds to setting the angle 
of E s to + (f 0 rather than - (p 0 (see Figures IV-4 and IV-7). The reason 
for the present approach to the adjustment problem is that for a weak 
station it is nearly impossible accurately to adjust the whip for a good 
null while the antenna is pointed in the reverse direction; this has been 
verified in the laboratory. 

Phase shift of the cardioid antenna system is plotted in Figures 
6-8 and 6-9, again using equations derived in Appendix IV. It should 
be noted that the Y-axis scales of the two figures are different. For 
0°< Ci< 180°, the phase shifts no more than 0.3 x 10“^ deg as the loop 
rotates. At a frequency of 20 kHz, 360° phase shift corresponds to 50 

O _ Q 

yUsec, so 0.3 x 10*"' deg would correspond to less than 10 ^usec; since 
normal phase variations of propagation produce changes of 1 or 2 jdsec 
on a phase plot of a received VLF station, these small changes would not 
be noticeable. 

An interesting aspect of Figure 6-9 is the phase change of -90° in 
the vicinity of (X = -90°. At 20 kHz, this would correspond to a phase 
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Figure 6-8. Cardioid Relative Phase Shift (0°<C(.<lS0 o ) 
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Figure 6-9. Caraioid Relative Phase Shift (-180°<<X<0°) 
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change of 12.5 ^tsec. It should also be noted that, while there may appear 
to be a discontinuity at Ot = 0° and (X = *±180°, due to the difference 
in T-axis scales on Figures 6-8 and 6-9 this is not in fact the case. 

Since it is proposed that a cardioid antenna system be used aboard 
ship in the VLF relative navigation system, there has been some concern 
that the sharp phase change in the vicinity of the cardioid null might 
cause difficulties. It is believed by the authors that this will cause 
no problems for the following reasons: (a) the relative broadness of the 
null, and (b) experimental evidence that the 90° phase shift in a typical 
cardioid system does not occur so abruptly as the theoretical curve of 
Figure 6-9 indicates. The latter will be discussed later in this sec- 
tion. The former is seen by examining Figure 6-7: if the loop is rota- 
ted 15° either side of the null, cardioid response is still down nearly 
35 db from that of the forward direction. For these reasons, it is felt 
that no special measures would need to be taken aboard ship to keep the 
antenna aligned with the direction of the transmitting station, other 
than a simple servo system driven by the ship’s gyro. 

Experimental Work 

The Tracor Model 611 Cardioid Unit was used at the Naval Postgraduate 
School. This unit has two inputs: one directly from a loop antenna, and 
the other from a whip antenna. Both antennas must be tuned as desired in 
the VLF band and matched to 50-ohm coaxial cable. The cardioid unit pro- 
vides adjustable phase shift and amplification so that the cardioid pat- 
tern may be obtained at any frequency between about 9 kHz and 30 kHz. 

The output of the unit is switchable between cardioid, loop only, and 
whip only. Although the cardioid unit is specified to be usable with 
only particular types of Tracor antennas, it was found to be usable with 
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an 18-foot whip antenna that happened to be on hand and a tuning network 
designed for 20 kHz. The complete cardioid system manufactured by Tra- 
cor comes with a remotely -tuned rotatable loop antenna and a whip 
antenna with fixed tuning. This system was used to obtain data pre- 
sented below. 

Radiation patterns for the cardioid system were obtained using WVL 
at 20 kHz. The length of a side of the square loop was 1.635 m, so data 
is directly comparable with the theoretical patterns of Figures 6-6, 6-7, 
6-8, and 6-9. The system was adjusted by setting whip amplitude and phase 
equal to that of the loop with the loop oriented toward Ft. Collins, 
Colorado (the location of WWVL) . It would have been preferable to set 
up a nearby low-powered (milli-watt or micro-watt) phase-stabilized 
transmitter to assure steady signal levels while the antenna was rotated, 
but this was not practicable. 

The amplitude data obtained is shown in Figure 6-10. Azimuth of 0° 
corresponds approximately to orientation toward the transmitting sta- 
tion. The numbers in parentheses below the abscissa correspond to angle 
of rotation 0L . Due to difficulty in reading the rotator unit, relative 
azimuth accuracy is probably *3°. Relative amplitude is accurate to 
better than il db. The two major sources of error here are variations 
in the received carrier level of WWVL and calibration of the receiver 
carrier level recorder output. 

It is interesting to compare Figure 6-10 with Figures 6-6 and 6-7. 
Considering that real antennas never perform quite so well as ideal ones, 
the similarity is striking. Overall, the null is somewhat broader than 
the theoretical; for example, the output is 35 db down at about 22° either 
side of the null. The 6 db beamwidth is 181°, which is the same as the 
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Figure 6-10. Cardiold Amplitude 



theoretical value, within the experimental error. The two points appear 
ing within the null at amplitudes of -32.5 db and -35.4 db indicate that 
the cardioid system was misad justed, with whip amplitude slightly great- 
er than loop amplitude. A careful examination of Figure IV-7 shows that 
this would also cause phase excursion to be slightly greater than A/4. 
This was verified experimentally, as seen in Figure 6-11. 

Figure 6-11 shows phase shift data obtained as the loop was rotated 
The T = 0 yUsec line is an arbitrary reference level. Phase data is 
probably accurate to 11 or 2 JJL sec, considering accuracy of the receiver 
and propagation changes that could have occurred while data was being 
taken. The smaller figure is probably acceptable, since data was taken 
during the daytime when the ionosphere was stable, over a period of 
about three hours. Angle of rotation ot is again indicated in parenthe- 
ses below the abscissa. The most notable difference between this curve 
and the theoretical data of Figures 6-8 and 6-9 is the broader width of 
the 90° ( A/4) phase change in the vicinity of the null, as noted ear- 
lier. The rubidium vapor frequency standard was used as the local phase 
reference for this data. 

There are various problems associated with using a cardioid antenna 
system in a shipboard environment. One is the problem of re-radiators 
located near the loop-whip system. This presents a problem in a receiv- 
ing system because of distortion of the antenna pattern, principally by 
resonant lengths of metal near the antenna. Unfortunately, it was not 
possible to investigate this problem at the Naval Postgraduate School. 

It is, however, the opinion of the authors that this problem would not 
seriously degrade the cardioid antenna’s performance in the VLF relative 
navigation system. 
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Figure 6-11, Cardiol d Pha 



Another problem is that of properly adjusting the cardioid unit. 

Whip phase and amplitude must be carefully matched to that of the loop 
at each frequency used, in order to obtain a true cardioid pattern. This 
may be difficult for distant stations whose long and short path signals 
are comparable. Careful attention should be given these adjustments, 
since depth of the null depends directly on how well the system is ad- 
justed. 
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7. Determination of Signal Strengths 

In this section values will be calculated for the relative strengths 
of signals received via the long and short great circle paths as a func- 
tion of the time of day and day of the year. This, along with the phase 
of both signals, is necessary in order to calculate the phase of the re- 
sultant signal. 

The relative strengths of the two signals can be found at any point 
on earth if the total attenuation on each path is known. The relative 
strength of the two signals is then just the difference of the short and 
long path attenuations. Upon reexamination of the electric field strength 
equation in Section 2 it will be seen that an assumption is being made 
here that all terms except the last two in this equation are identical 
for both propagation paths. It seems only reasonable to assume that most 
VLF stations have a symmetrical radiation pattern and that the amount of 
excitation in two opposite directions will be equal. 

Attenuation rates at VLF, mentioned in Section 2, vary from about 
one to five db/Mm depending on the time of day and propagation path con- 
ditions. References 5, 8, 9, 18, and 28 contain many of the curves of 
attenuation rates versus frequency which have been derived from experi- 
mental data. These curves will be used extensively here. 

The procedure used to obtain expected attenuation rates at any par- 
ticular frequency consists simply of sifting through the data present in 
these references and extracting the attenuation rates given for different 
conditions at this frequency. From this information estimates of the 
attenuation rates for all possible path conditions and times of day can 
be made. 

As an example, at 13.6 kHz the following values are given in the 



63 



references stated. 



OC (db/Mm) 


conditions 


reference 


3.6 


land, day, east to west 


Taylor [ 9 ] pg. 1933 


2.8 


sea, day, east to west 


Taylor [9] pg. 1933 


2.6 


day 


Chapman and Macario C$3 
Pg. 930 


2.0 


sea, night, east to west 


Taylor and Lange £83 
pg. 609 


1.8 


day, sea, west to east 


Taylor [93 pg. 19 33 


1.6 


night 


Chapman and Macario [5^ 
pg. 930 


3.5 


sea, day, east to west 


Watt and Croghan £18^ 


2. 35 


sea, day, north to south 


Watt and Croghan [18^ 


1.95 


sea, day, west to east 


Watt and Croghan [18^ 


1.3 


night, sea 


Watt and Croghan £l83 


1.0 


sea, night, west to east, 
15-20 kHz 


Davies £28^ pg. 425 


= 1.0 to 1.5 


sea to land change, 
15-20 kHz 


Davies Q2 pg. 425 


A «= 1-0 


night to day change, 
15-20 kHz 


Davies [28j pg. 425 


Ac <= l.o 


W-E to E-W change, 
15-20 kHz 


Davies [28^ pg. 425 


From these attenuation rates (and changes in attenuation rates) and 


other discussions of 


attenuation rates such as those by Crombie CO and 


Westfall [2l] reasonable estimates of attenuation 


rates for 13.6 kHz can 


be made. These will 


be in the form of one basic 


attenuation rate plus 



additional factors which depend on the direction of propagation, amount 
of daylight on the path and amount of land on the path. 

The values chosen are: 
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0^ = 1.0 db/Mm 

AOT n -d = 1 -° d bAfa 

AO<s_i = 0,8 db / Mm 

Ao<d N _s = 0,4 db / Mn 

^°*dE-W = 1 - 55 db/Mm 
The change in attenuation 
approximated by: 



basic rate (night, sea, W-E) 

night to day change 

sea to land change (average land) 

W-E to N-S direction change 

W-E to E-W direction change 

e as a function of direction can be 




where p is the direction of propagation with respect to 090° true 
(0£^^180 o ). For north to south propagation 0 = 90° and for east 
to west propagation 0 = 180°. 

The attenuation rate over any path can be calculated from the fol- 
lowing equation: 



* = <*b + P D ^*n-d + P L*<*s-l 
where P D = fraction of the path which is in daylight 
P^ = fraction of the path which is over land. 

A typical plot of attenuation for a 24 hour period is shown in 
Figure 7-1. The path here is Trinidad to Monterey (short path); 0 = 
145°, P L = 0.6 and P D is obtained from the computer program described in 
Appendix II. The total attenuation is found by multiplying 0< by the 
path length obtained from the geodesic program (Appendix I). The cal- 
culated attenuation rates vary from 2.48 db/Mm for the entire path in 
darkness to 3.48 db/Mm for the all-daylight path. 

Figure 7-2 shows a comparison of the long and short path attenua- 
tion from Trinidad to Monterey. The long path here is considered to be 
in the opposite direction at the receiver and transmitter from the short 
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Max. =23. 37 db 
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Figure 7-2. Calculated attenuation of 13.6 kHz signal on 
Trinidad to Monterey long and short paths for 
1 Oct. 1966 

Short path length = 6658.85 km 
Long path length = 33391 • ^8 km 
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path direction of propagation. 

It should be noted that the attenuation rates used here for land 
paths are for average land. Large areas of very high or very low conduc- 
tivity land would require alteration of these values. Similar care must 
be exercised for propagation over polar ice caps, where the conductivity 
is lower than that of sea water. 

One factor not mentioned previously is the loss due to discontin- 
uities of the earth conductivity along the path (L^ n in the electric 
field strength equation). Westfall C 21 D mentions this factor briefly 
as having a magnitude of three to six db at times. Because of the lack 
of information on this factor it has not been taken into account here. 

The method shown in this section could be used at any other fre- 
quency of interest in the VLF range to obtain similar attenuation rates. 
Some complicating factors may be encountered, however, which could make 
the calculations more difficult. For example, at 16 kHz it was found that 
the night to day difference over land is greater than the night to day 
difference over sea. 

Comparison of Experimental and Calculated Values 

Figures 7-3 and 7-4 show two more curves of calculated attenuation 
for the Forest port and Haiku Omega stations and Figures 7-5 through 7-7 
show the observed relative carrier levels of these three stations as ob- 
tained from a Tracor Model 599G VLF phase tracking receiver. No measure- 
ments were obtained on any long path signals or on the Aldra short path 
signal because the carrier levels were below the level of control of the 
receiver AGC, giving a reading of zero on the relative carrier level out- 
put. 

Since no attempt was made to calculate the signal strength of the 
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UT 

Figure,. 7-3* Calculated attenuation of 13*6 kHz signal on 

Forestport to Konterey short path for 1 Oct. 1966 
Path length = 4011.74 km 
Min. = 12.90 db 
Max. = 16.91 db 




UT 

Figure 7-4. Calculated attenuation of 13.6 kHz signal on 
Haiku to Konterey short path for 1 Oct. 1 966 
Path length = 3856.50 km 
Kin. = 4.16 db 
Max. = 8.02 db 
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UT 

Figure 7-5. Observed relative carrier level of Trinidad 

13.6 kHz at Monterey, 30 Aug. to 2 Sept. 19 66 

10 db/l- in. vertical scale 

night average - day average = 4.04 db 




Figure 7-6. Observed relative carrier level of Forestport 
13.6 kHz signal at Monterey, 5 Sept, to 8 Sept- 
1966. 

10 db/ln. vertical scale 

night average - day average = 2.47 db 




Figure 7“7» Observed relative carrier level of Haiku 13-6 

kHz signal at Monterey, 2 Sept, to 5 Sept. 1966 

10 db/jjr in. vertical scale 

night average - day average — 0.3 db 
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various stations and this information is not readily obtainable from the 

Tracor receivers, the only quantities that can be compared meaningfully 

are the difference between the carrier levels for the entire path in 

darkness and in daylight and the day-night attenuation difference. While 

the attenuation decreases at night and the carrier level increases during 

the same period the difference between the day and night values should 

be equal since they are complementary quantities. (This is easily seen 

in the electric field strength equation of Section 2 where C<d/10^ is 

the attenuation and is the carrier level.) 

z , n 

A comparison of the calculated and observed values of the night -day 
difference is shown below. 



Station 


Calculated difference 


Observed 


difference 








(a) 


(b) 


Tr in id ad 


6.65 db 




4.04 


6.4 db 


Forestport 


4.01 db 




2.37 


7.5 db 


Haiku 


3.86 db 




0.3 


4.5 db 


The values 


listed under (a) were 


observed 


with the 


receiver servo 


constant set at 150 seconds and those under (b) with a 15 second 


constant . 


The lower differences 


observed 


with the 


150 second time 



constant, corresponding to Figures 7-5 and 7-7 are apparently because 
of the erratic nature of nighttime phase. 

Calibration of the signal level scale used on the recorder was very 
difficult, leading to some possible inaccuracies in these measurements. 
(Calibration was complicated by the short duty cycle of the transmitting 
stations; 10%. ) Thus while these measurements do provide some measure of 
confidence in the calculated attenuation rates, it is believed that more 
comprehensive measurements of the carrier levels of these signals would 
be necessary if an accurate check of these attenuation rates was desired. 
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8. Calculation of VLF Phase Velocities 

Figures 2-8 through 2-11 show several curves of calculated phase 
velocities for VLF signals for different conditions which could be used 
directly to obtain the desired phase velocities. However, it would be 
more convenient and accurate to make use of an equation for phase velocity 
which takes into account the frequency and conditions of the ionosphere 
and the earth. Such an equation is available in the Omega System book 
[3l], page 3-20: 




0. 36 




(2 Ifn - 



w - yO 



V o ]2 

4JTV2^fh i J 



Here Vp = phase velocity of the VLF signal 
V Q = free space velocity of the wave 
= 3 x 10^ m/sec 

0.36 = a constant chosen to give a good fit in the 10-20 kHz 
frequency range 

h^ = ionospheric height for which the conductivity (tT) = 1CT 6 
mho s/m 

a = radius of the earth = 6378.155 km 
g = phase shift upon reflection at the ground 
= phase shift upon reflection at the ionosphere 
n = mode number (= 1 for primary mode) 
f = frequency 

Some values are given for the phase shifts that can be expected; 

= -160° for daytime propagation, -165° at night, -165° for daytime 
propagation to the east, and -155° for daytime propagation to the west. 
Assuming that the -160° phase shift is for north-south propagation, which 
seems reasonable, the values for nighttime propagation to the east and 
west would then be -170° and -160° respectively. The ground phase shift 
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is zero for sea water and the values given for land are -0.09 radians 
during the day and -0.13 radians at night, for a ground conductivity of 
10 ~^ mhos/m. 

Using the given ionospheric heights of 68 km during the day and 86 
km at night (derived from an exponential conductivity profile and taking 
the height for which KT = 10“^ mhos/m) the phase velocity curves of Fig- 
ure 8-1 were calculated. A comparison of these curves with similar 
curves in Section 2 (Figures 2-9, 2-10, 2-11) shows a very close corres- 
pondence. Several experimental phase velocities from Figure 2-14 are 
also plotted here for further comparison. 

The phase velocity at any point along a path is: 



Vp (x) 



2 TT f 
d Q /dx 



where & is the phase in radians at the point x on the path. The total 
phase along the path is then: 

e - 

If the phase velocity is constant this becomes: 



_ 2 TTfd 
V P 

Using a method similar to that employed in determining the path 
attenuations of Section 7, an average phase velocity can be obtained 
for each path. It can be calculated by the following equation: 



'pi, * <*> 



+ P n ^V „ „ + P„AV , 

pN-D S pi— s 



where 



A V, 



pd 



V . = basic phase velocity; phase velocity for night, land, 

W-E from Figure 8-1. 

(0) = increase in phase velocity due to direction of propagation. 
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P D = proportion of path in daylight. 

AVpN D = increase in phase velocity for daylight propagation. 

Pg = proportion of path over sea. 
av pi- g = increase in phase velocity for propagation over sea. 

The increase in phase velocity with direction can be obtained from 
the following equation: 



where 0 is the direction of propagation, measured as in Section 7. From 
Figure 8-1, the differences between the phase velocities for E-W and W-E 
propagation are different for daytime and nighttime propagation. * y pd 
for daytime can be approximated by the average daytime differences in 
phase velocity with direction and the nighttime value can be obtained 
similarly. This can be expressed by the following equations: 

^ V pd = ^ (V psDEW ~ V psDWE + V plDEW " ^lDWE 5 da y time 
^ V pd = ^ (V psNEW " V psNWE + V plNEW ~ V plNWE^ nighttltne 



Here V p S p)EW = P hase velocity for sea, day, E-W 

'VplNWE = phase velocity for land, night, W-E and so on. 

Using reasoning similar to the above it is found that two values 
for D must be used (one for E-W propagation and one for W-E) and 

one value for can be used. These are given by the following 

equat ions : 

^ V pN-D = ^ (V psDWE " V psNWE + V plDWE " V plNWE ) ’ W_E 

^ V pN-D = ^ (V psDEW ' V psNEW + V plDEW " V plNEW ) » E_W 

^ V pl-s = ^ (V psDWE " V plDWE + V psDEW " V plDEW 
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+ V psNWE " V plNWE + V psNEW " ^lNEW^ 

The range of the phase velocity changes being averaged to obtain 
A v pi-s a ^ out ^ x 10”^ V p /V Q . The range of phase velocity changes 
being averaged for * V pN- _ D and AV pd i$ smaller, indicating that the 
error in AVp^_ s is the largest that would be encountered. 

The curve shown in Figure 8-2 is a typical 24 hour plot of phase 
delay calculated by the above method. This is for the Trinidad to Mon- 
terey short path with the same parameters as used in Section 7(0 = 
145°, P g = 0.4, P D from computer program of Appendix II). The 13.6 kHz 
phase velocities used to calculate this curve are listed below. 



V psEWE = 
V psDEW = 
V psNWE 
V psNEW = 
V plDWE 
V plDEW = 
V plNWE = 
V plNEW = 



1. 00002 13V 0 
1. 000427 3V 0 
0.9974367V 

o 

0.9976842V 

o 

0.9998199V 

o 

1.0002154V o 

0.9972606V 

o 

0.9974986V q 



Comparison of Experimental and Calculated Phase Delays 
Figures 8-3 through 8-6 show both calculated and observed 24 hour 
phase delay curves. The average of the observed daily curves is shown 
shifted by 1.5 inches (15 or 30/ksec) for convenience and ease of com- 
parison with the calculated curves. The total phase delay of the signal 
from transmitter to receiver is not easily obtained experimentally; thus 
only the day-to-day changes in the phase delay are compared here. 

Figures 8-5 and 8-6 show a very close correspondence between calcu- 



76 



lated and observed curves. It is interesting to note that there appears 
to be a time lag in the observed phase delay during periods of rapid 
change, which could be due to the long time constant used in the receiver 
(150 sec). The lag also seems more pronounced on stations with lower 
signal strengths at the receiver. 

Figure 8-4 shows differences in the calculated and observed curves 
which seem to point towards some error in the calculations. The ampli- 
tude and shape of the observed diurnal shift are believed to be approx- 
imately correct. 

Figure 8-3 is a good example of a diurnal shift that is not a trap- 
ezium. Propagation over a polar path, as is the case here, yields 
diurnal shifts similar to those shown in Figure 8-7, and are markedly 
different during different times of the year. The diurnal shift on polar 
paths is discussed more thoroughly by Blackband Cl2j. In Figure 8-3 the 
poor correlation of the calculated and observed curves may be partially 
due to the weak signal strength of the Aldra station at Monterey. The 
phase-tracking receiver did not seem to "lock on" to the received phase 
until a large change in phase delay occurred, such as encountered during 
sunset along the path. This can be seen during the initial portion of 
the readings, started at 1900 on 23 September. 

Figure 8-8 shows a comparison of an experimental phase velocity 
curve, obtained by Thompson et_ al [14] in Sydney, Australia at 16 kHz, 
with a calculated curve. The data used for the calculated curve is 
approximate (position of the receiver, proportion of land on the path, 
etc.) and in view of this the close correspondence of these two curves 
is remarkable. 
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Figure 8-2* Calculated phase delay on Trinidad to Monterey 
short path at 13.6 kHz- 
Minimum value is 2219.5 • 759 /^-sec delay 
10/u.sec/inch vertical scale 
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Figure $-3'. Aldra to Monterey phase delay, short path 
13.6 kHz, 23-26 September, 1966 - 
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Klgure' 8^-4. Haiku to Monterey phase delay, short path 
13*6 kHz,. 25 July- 1 August, 1§66, 
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Figure 8 - 5 . Forestport to Monterey phase delay, 
13.6 kHz, 6-9 October, 1966 



short path 
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Figure: 8-6. Trinidad to Monterey phase delay, short path 
13*6 kHz, 11-14 October, 1966 
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Figure 8-7. Typical diurnal shifts on a polar path. 
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Figure 8-8. Observed and calculated phase of GBR at Sydney, 
Australia for October 1961 and October 1966 
respectively. 
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9. Calculation of Resultant Signal Phase 

As discussed in Section 3, the amplitude and phase of a signal 
consisting of long and short path components can be obtained by the 
addition of two phasors. In the navigation system proposed in this 
paper the phase of the resultant signal at a receiver is of primary 
interest. The amplitude of the signal is of interest only insofar as 
it is sufficient to allow phase tracking. 

The phase of the resultant signal can be found from the phase of 
the long and short path signals and their relative amplitudes or atten- 
uations. First take the magnitude of the phasor in Figure 3-3 to be 

AX 

unity, or zero db. Then the magnitude of E 2 is the relative attenuation 
of the short path signal with respect to the long path signal. It can 
reasonably be assumed that when ) E 2 I ^ 0.1 the phase of is not appre- 
ciably affected by E^ (for = Q « 1 ^ 90° . = 5.7° in Figure 3-3). 

The relative attenuation for this case is given by the following 
equat ion: 



A = 20 log 



M 

|E1| 



20 db. 



Thus the long path signal must be taken into account whenever the 
strengths of the two signals are within approximately 20 db of each 
other. 

To obtain the phase of the resultant signal, the phase delay in 
microseconds, calculated in Section 8, must be converted to an angle. 
This is easily done by converting to wavelengths, where one wavelength 
equals 1/f seconds (at 20 kHz, 1/f = 50/Hsec). One wavelength corres- 
ponds to 360° and 0 is then the difference between the phase delays on 
the long and short paths. The angle 0 is then given by the following 
equation : 
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tan 0 = 



| E 2 1 sin 0 
1 + | E 2 | cos 0 

This angle could also be found by adding graphically the phasors and 

e 2 . 

It must be remembered here that the total phase delay of the re- 
sultant signal is equal to 0 plus the phase delay of the short path sig- 
nal . 

Figures 9-1 through 9-5 show the calculated curves of long and short 
path attenuation, and short path and resultant phase at Sydney, Australia, 
and Ellington, New Zealand. The transmitter station used was GBR at 16 
kHz and the day was assumed to be 2 3 September 1966. 

In Figure 9-2 the largest difference between the short path phase 
and the resultant phase is slightly less than one microsecond and the 
difference in attenuations at that time is 20 db. This verifies the 
earlier assumption that the resultant phase differs from the short path 
phase (by an appreciable amount) only when the strength of the long path 
signal is within 20 db of the short path signal. 

Figure 9-4 shows a case where the resultant phase is considerably 
different from the short path phase. The largest difference here is 12.1 
/Asec (at 0300) or about 70 electrical degrees, and does not occur at the 
time when the long path signal is strongest (long path signal is 0.6 db 
weaker at 0 300, 1.8 db stronger at 0500). From about 0800 to 1800 there 
is approximately 20 db of separation between the two signals and the re- 
sultant phase is very close to the short path phase as expected. 

In Figure 9-5 the latitude of Wellington was changed by four min- 
utes, lengthening the short path by about seven km and shortening the 
long path by the same amount. The total phase delay on the short path 
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Figure' 3 -i> Attenuation on GBR to Sydney short and long paths 
at 1 6 kHz. 
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Figure 9-2. Short path and resultant phase of GBR at Sydney 
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Figure 9-3* Attenuation on GBR to Wellington short and long 
paths at 16 kHz. 
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Figure 9-&. Short path and resultant phase of GBR at Wellington 
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Figure 9-5*. Short path and resultant phase of GBR at 
Wellington. (Short path 7 km longer) 
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is increased slightly as a result of this, but the diurnal shift remains 
the same. The resultant phase, however, is considerably different, show- 
ing that a change in position has a larger effect on the resultant phase 
than on either the long or short path phase alone. 

No experimental data has been found with which these calculated 
curves could be directly compared; however, Thompson et_ al . [14] and 
Crombie |[7] have published articles which discuss this phenomenon and 
include some experimental evidence taken at Sydney and Wellington. 

It should be noted that the resultant phase curves as shown here 
could not be observed by using a loop antenna. The 180° phase reversal 
with a reversal of the direction of the incoming signal would cause 
some different resultant phase to appear. This could be taken into 
account in the calculated curves, however, by adding 180° to the phase 
of the long path signal before adding the phasors. 

The concept of phasors as discussed here provides another means of 
describing the phenomenon of cycle-slipping. Figure 9-6 shows the am- 
plitude and phase of some station over a 24-hour period as represented 




Figure 9-6. Phasor representation of a 24-hour 
phase and amplitude curve 
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This figure is drawn assuming a trapezoidal diurnal pattern of both 
phase delay and amplitude, both being greater at night. Thus % repre- 
sents the daytime phase and amplitude, B represents the nighttime phase 
and amplitude, is the diurnal phase shift, and |B| - |A| is the night- 
time increase in amplitude. 

Cycle-slipping has occurred if the curve traced by the phasor en- 
closes the origin, as shown by the dotted line in Figure 9-6. This 
could be due to long path interference, loss of signal (in which case 
the magnitude of the vector is zero), or several other causes. The 
case for which the diurnal shift is greater than one cycle ( 2 TT ) should 
not be confused with cycle-slipping as shown here. Figure 9-7 shows one 
possible curve of this type and it can be seen here how easily cycle- 




A rapid change in phase, as sometimes experienced during sunset or 
sunrise, could cause the carrier level in the receiver to be greatly re- 
duced for a time, particularly if a long servo time constant is used. 
This could result in a phasor trace similar to the one shown in Figure 
9-8, with an increased probability of cycle-slipping. This rapid phase 



93 



change and resultant cycle-slipping is typical of north-south propagation 



paths . 




Figure 9-8. Phasor representation of loss of carrier 
level due to rapid phase change 
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10. Conclusions 



Attenuation rates and relative signal strengths of long and short 
path signals can be determined easily. Ihe accuracy of the method used 
here has not been determined, however, because of a lack of dependable 
experimental data. 

It has also been demonstrated that the phase of a signal can be 
predicted for a typical 24 hour period, both with and without long path 
interference. Experimental data has shown the short path phase calcula- 
tions to be very accurate for some cases when three-to five-day averages 
of observed phase are considered. Lack of experimental data prevented a 
check of the calculated phase in the presence of long path interference; 
however, it was shown that when this type of signal is present the phase 
is highly sensitive to a change in receiver position. 

Additional experimental data, possibly from some location other 
than Monterey, California, would be highly desirable and would probably 
lead to some modifications of the parameters used in the calculations 
presented here. It should be emphasized that these calculations are a 
first iteration in an attempt" to predict, at any point on earth, the 
phase of a VLF signal. 

A cardioid antenna system should be suitable for use aboard ship 
in a VLF relative navigation system. A simple servo system driven by 
the ship's gyro, and provided at intervals with the true bearing to the 
desired VLF station, should be adequate to keep the cardioid pattern 
oriented with sufficient accuracy. With proper adjustment of the car- 
dioid system, interference from the long path signal may be eliminated. 
In case of a location at which the long path signal is always equal to 
or greater in strength than that of the short path, however, it may be 
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impossible to adjust the cardioid system properly, and hence it could 
not provide the desired amount of long path signal attenuation. 
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APPENDIX I 



GEODESICS AND LONG PATHS 

A. Geodesics 

The great circle distance between two points on a sphere is easily 
calculated; however, determination of the length of a geodesic (analo- 
gous to a great circle) on an ellipsoid is more difficult. 

Kirkland f2C0 states very concisely the simplest method of obtaining 
the length of a geodesic. 

Given two points ( Ap, (0 2 * ^ 2 ^ their respective para- 

metric latitudes, ^ and yg 2 as shown in Figure 1-1. 

0 = geodetic latitude 
/ 3 = parametric latitude 
= equatorial radius 
b = polar radius 



dz 

<*y 



tan 0 - - 

z = a cos^ 

y = b sin^ 

dz = - a tan £ 
dy b H 

tan ^ tan (fi = c 

sec^ = 1 + c^ tan^ 0 

COS = 1 




sec 



7T 



cos 



Ycos 2 



£ 



+ c^ sin^ 



sin^ = Vl - cos ^ 1 



c sin 






^jcos^-0 + c ^ sin^jb 1 
The arc distance between the two points (x) is then found from: 
cos x = sin (£ y sin 2 + cos ^ ^ cos y$ 2 cos ( A 1 ■ A 2 ^ 
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and the great circle distance is D = a x . 

Now using the correction factors derived by Lambert Cl3 the length 
of the geodesic is: 

d = D - A q P - B q Q 



where 



(a - b) (x - sin x) 
4(1 + cos x) 



Q = 



(a - b) (x + sin x) 
4(1 - cos x) 



A q = (sin (Z y + sin 

B q = (sin /3 \ - sin @ 2 ) 2 

The accuracy of this method is, according to Lambert, better than 
i 0.05 km but breaks down near the antipodes (i.e. for x^>175°). This 
method was checked against the examples in an article by Rainsford CO 
with the largest error being about 0.85 km. (It should be noted that 
Rainsford specified the position of the points to within 0.0001 second.) 
B. Long Path Distances 

The method used in this paper to obtain the long path distance be- 
tween two points is to find the length of the perimeter of the ellipse 
on which the two points lie (analogous to the great circle on a sphere) 
and subtract the geodesic length from this distance. 

The major axis of this ellipse is the equatorial radius of the 
ellipsoid and the minor axis can be found as follows: 

@ = parametric latitude 
x = arc distance between the two points 

cos (X = cos A - cos X cos C 
sin x sin C 
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Pole 



sin /3 2 ~ cos x s ^- n /3 l 
sin x cos ^ 




Eau^roR 

Figure 1-2. Ellipse of geodesic 



v s> /o fsin/3o - cos x sin^n 

cos o = sm (X cos Q = cos p -\/ 1 - 1—4 1— i—L\ 

^ 1 r 1 \l [ sin x cos ^3 l J 



T 



Here gamma is the parametric latitude of the northernmost point on the 




The minor axis of the ellipse of interest is then: 



- 4^7 



= -\/b 2 + (a 2 - b 2 ) cos 2 tf 



The length of the perimeter of the ellipse described by a and r is 
found by evaluating the elliptic integral: 



E(k) 



jVr 



2 o ' 
k sin x dx 



where 



* V 1 - 
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The perimeter is then: 

P = 4 a E (k) 

and the long path distance is: 

a L = p - a 

where d is the geodesic length. 

The FORTRAN 60 program used to compute these distances is shown on 
the following pages along with sample results. 
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PROGRAM GEODESIC 



I 







LU 




















1 








1 
















j 


i 

i 

i 
















t 






X 




























i 




























! 


1 










►H 


























































O 










X 






* 
























X 




























•— * 










cr. 






o 
























X 


























i 


CO 




f 






UJ 




o 


2 
























< 




























X 




1 


LU 




x 


LU 




*— * 


















CO 






X 


2 


























3 






X 




LU 


X 


o 


2 


















X 








o: 






















i 


o 




i 


X 




O 


X 


X 


LU 


















CO 


2 




X 


•— * 


























X 








LU 






LU 


X 




















— • 




3 


1 — 


























CD 




1 


2 


o 


CJ 


CO 


X 


X 














i 


2 


3 






o 


•— • 
























i 


I 






O 


2 




•— • 


X 


<C 




















O 




X 


00 


































< 


o 




CO 


3 


















x 






•— < 


O 
































CO 


f— 


2 


00 




X 
















•H 


3 


CO 




o 


X 


1 






























1— 


CO 


<< 




LU 
















! 




X 


LU 






j 


! 






























2 








X 


LU 
















CO 


X 


O 




) — 


X 


1 














* 


















3 


•s 


CO 


i — 


X 
















LU 


3 


O 




< 


X 


[ 


+ 










| 


















o 




CO 


3 




X 
















cr 


O 


X 




X 


X 




LU 




























CL 


1 — 


•s 




•> 


















LU 




o 




X 






> 






























CO 


< 


3 




X 
















X 


3 




i 


CD 


X 




•—4 




























O 


LU 




C 


c 


o 
















2 


3 


X 






O 


F 

i 


X 




























3: 


r- 


• 


X 


\ 




X 
















X 


X 




X 






x-« 




























i— 


cr 


l— 




CO -J 


















2 


X 




X 


X 


i 


CO 






























o 


x 


-J 




<: 
















X 


•— * 






X 


3 


O 




























2 : 


X 


LU 


< 




CJ 


i n 


00 








2 




o 


X 


3 


1 




3 


X 




























LU 


CO 


CD 




o 


o 


in 


CM 


*: 












X 


2 




X 


x 


X 






























LU 




X 


X 




X 


• iH 


tH 












CO 




< 




o 


►— * 


o 


CO 
























X 




■v 




< 


o 


X 


CL 


x rv 


00 


CO 




in 


ro 


j 


I - 


• 




| 




CD 




•— 1 
























3 




►— 


CJ 3 


X 






\ X 


OJ 


oo 




in 


rv 




2 


o 


(A 








< 


























o 




LU 


>— « 




<. 


LU 


O 


ro 


LA 


ro 




x 


rv 




— • 


X 


2 




2 


o 


3 


• 
























X 


J 


X 


CO 


• 


3 


X 


LU 


. 


T“f 


* 


o 


• 


• 


ro 


o 


< 


< 


j 


o 


3 


O CD 




















! 


i 


— • 






LU 




O 


X 


X 


i hvi • 


00 


• 


CO 


NO 




X 


o 


— ■ 






X 


2 
























o 




x 


3 


* 


LU 








On 




rv 




in 


CO 






O 




X 


X 




o 




























o 


o 


3: 




X 


X 


x ro 


On 


PO 


o 


ro 


ro 


o o 


< 


< 




< 


X X 


3 








i 


J 














X 




2 


UJ 




LU 


o 


X 


NO 


CM no 


CM 


nO 


o 


CM 


2 


X 


X 




2 X 


x 












| 














< 


j 


< 


o 


>- 


X 




1 J— 


ii 












<X 


< 








: 


X 


X 








I 
















X 




X 




CD 


I — 


o 


it 


II 


ii 




1 11 


ii 


ll 


1 


a 


2 




3 Q 




CO 










i 












i 


X 




CO 


LU 








•>> 


o 














X 




i 




CJ 


2 


CO 


X 










[ 














CD 






X 


Q 


CO 


X 


X 


<c 


X 


c 


X 


a 

V 


X 


X 


UJ 






X 


< 


X 


3: 










i 












i 






3 


X 


LU 




< 




LU 


1 












X 


2 


X 


1 


Q 




00 












i 












X 










X 




X 


Q 


CO 
















O 


2 






X 


CO 


•s 






















X 






LU 


CD 


O 






2 


<c 


i 


o 




o 






X 








X 


o 


O X 






















< 






X 


7 


— 1 




UJ 


<i 


o 


j 


o 




NO 






o 


2 


O 




X 


X 


rr 


•> 






















CL 






Y— 


o 


LU 


► 


X 








►— « 




ON 








Q 


X 




X 


< 


CJ 


•— * 


/ 














1 










J 




3 


r> 


o 


X 


< 


X 




X 




H 






• 










3 




i — , 








\ 












1 






to 


< 


LU 


LU 






o 




X 










CJ 




X 


i • 


CO 


CJ x 










! 


i 


















uj 




o 


’A 


CO 


V 


— • 




X 










<: 


• 


CJ 


2 




X 


3 


CO 










i 

i 


\T 


nO 








j 






b— 


X 




3 


•— « 


QD 


X 




X 




Q 






o 


CJ 


< 






(J 


o 










O 


I 


1 










X 




3 


*— 


CO 








3T O 




CO 




O 




i 


1 I 


LU 


X 




X 


X 


X 


X 










X X 


IU 








1 


X 




X 




2 


LU 


o 


O 


o 








— • 








CO 




a 


3 X 


' • 












in 


NO 


X 








j 


< 




X 


<r 


O 


CD 




X 


2 




3 




X 






Q 




X 


2 


X 


X 


CJ 


CO 


o 








CM 


NO 


X 










X 




o 


LU 






o 




— ' cc 


<r 


i 


X 






2: 


3 


o 


C 


X 




! 




ro, 








On 


00 


oo 














o 




J — 


O 


2 


X 


LU 




2. 


1 


X 






< 


2 




1 


3 


o 


X 




in 








CM 


CD 


NO 










CD 






UJ 


<t 


i — 


<j: 




- X 




O 




X 








< 


• 


CO 


O 


X 


< 


X 


NO 






ro 


Od 


to 








2 




X 


X 


X 






CJ 


~ X 




•— * 


i 


00 






| • 




CD 


X 




X 


3 


cm 








in 


CO 


X 










O 




<c 


1 — 


3 


O' 


CO 


X 


CD CO 




X 










X 


• 


2 


3 


X 


X 


X 


3 


o 


in 


ro 


1 


-r 


00 


CC 








3 




CSC 




LU 


o 


3 


X 






< 




3 






< 


2 


O 




3 


O CD 


X 


in 


in 








CD 


M* 








X 






5 


X 




— « 


— • 


CO 


1 3 




2 


i 

; 


X 






3 


•— * 


3 2 


O X 






tH 


v - 1 r-v 




i 


o 


CO 








X 






o 


O 


LU 


X 


O 




LU 




X 




O 








2 






X 


XXX 


m- 


• 


• 




CD 


• 


• 








On 






x 




X 


X 


< 


2 


C CO 




UJ 




X 






CO 




a 


3 




CO 


X 


< 


X 


co 


NO 




• 


CM 


M' 


3 






O 






X 


LU 


X 


X 


X. 


,x 


^ CO 




X 




X 






•— * 


CO 


2 


< 


O XX 


3 






in 














O 


'w' 








3 




X 




lx 


\ LU 




2 




o 








LU 


< 


CJ 




X 






ro 


ro 


ro 




II 


II 


II 


2 


X 




X 




f 


CO 


< 


CD 


CO 


X 


X 


c x 




•— * 




00 






X 


LU 






X 


1 


X X 




NO >0 


<r 








O 


CO 


X 


< 




x> 


•— . 


X 


2 




< 


o 
















3 


X 


• 


X 


c 


X 


X 


X 


n 








3 


2 


(A 


2 


a 


2 


2 




2 


X 


or 


•— » 


LU 


3 




II 












X 


CD 


X 


3 


X 


X X 


X 


1 


II 


ii 


CO 


3 


3 


3 


2 


—+ 


*— • 


X 




*— 


X 


X 


CO 


X 


o 


CO 
















2 


LU 


<f 


< 


X 




, X 


O 








ii 






<C 


O 


3 X 


o 




X 




CO 


3 


X 


X 


; 


X 














— • 


Q 


3 


2 


CJ 




3: 


2 


X 


c 


CO 


O X 


X 


X o 


X X 


X 




X 






































j 




























ro 




























i 




































1 

| 


NO 






o 


j 

o 


CJ 


o 


o 


o 


o O 


o 


o 


o 


o 


o 


<_> 


O 


O 


O 


O 


! 

j 

CJ o 


CJ 


O 














j 




X 


} 

| 






i 



I 



I 



104 



3 

<z 













; 






































1 




























































i 


1 

i 




V— 






















































i 


i 




CO 






















































1 






3 










o 


















































X 










CM 




























































-i 


















































•k 










xj 


















































»<-» 




j 










CO 














































CO 










CO 




<c 














































X 










CM 




M- 


















































I 






~ 


















































•k. 




i 






X 




•k 














































UJ 




1 






X 




<=> 














































Q 




j 






•k 




• 








, 






































3 










X 




ro 














































►— 










CM 




X 
























































■— « 




CM 














































h- ' 










X 




« 














































< 










OL 




o 














































3 










o 




in 


CO 




CO 






















3 


















O 










CM 




X 


tH 




CM 






















CO 








y— * 








— * 




j 


• 




— « 






•—4 




•—* 










1 












o 
















X 






o 




X 






X 




x : 






















o 










CM 


#— kk 






V-' 






2 




CL 




o 


X 




Xj 






















*1 


! 










CM 




X 






o 








• 


* 




* 




















*»4k 


1 








CO 








UJ 




i 


3 




CO 




ro 


CO 


CO 


CO CO 




















CM 






<— V 


z 


00 




o 




1 






xH 




X 


Q 


tH 


Gi 


CM 










• 










W| 


j 








*— t 


z 




o 










3 




CM 


<t 


-J 


c 


3 




















CD 








X 


CO 








UJ 




I 






fX 






X 


rxl 


cr 


r-M 




















co; 






X CO 




(O 




3 




— «k 






•k 




•kj 




4 


1 


4 






• 














O 




CO o 


4 


1 
















X 




o 


4 


CO 


+ 


CO 




















O 






o o 




1 






UJ 


*** 




• 


•*~<k 


tH 


o 


•I 




Q 


! 


Q 




















*J 


! 




a 


1 










X 


tH 


CM 


3 


CM 


3 


CO 


in 


X 


< 


X 


< 








1 














i 






! i 


tH 


<— >k 






h- 




w 


u 


s-/ 


NOIl 


H 


X 


CM 


X 




















t-l 


, / 




4 




•W* 


tH 


i 






4— H 




UJ 


CD 


* 


X 




•— «* 




•—4 






















3 


l 

i 






i • 


CD 








CO 


X 


X 


Q 


c/r o 




o 


X 


+ 


x ; 


4 




















00 




« 


tH 




CD 




X 


a. 


X 




o 


tH 




»; 


X 




x5 






















CO 






tH 


1 


►— ♦ 


2 




o 


w 


k— ^ 




o 


j o ro 


4 


X 


*1 X 










CM 










O X 




w 


; 4 


CO 


•— * 


> 




t-H 


X 


X 




•» 


M CO 


X 


X 


tH 


X 


CM 










X 










O CO 




4 


• 


w 


CO o 




h- 


X 


X 




3 


•> 


o 


CM 


o 


3 


o 


3 










a 












O 






M* 


4 




¥ 




< 


<C 


< 




CM 


CO 


X 




<c 


fsi 


<t 


fSI 










z 










+ 


o 


/ 


M" 






4 


o 




h-‘ 


) — 


' — 




^ t-e 




CD 


X 


4 


X 


4 










CO 










1 


4 


/ 1 




X 






CD 




3 


4 


4 




CD 


t — » 


•* 


• 




X 


i 


>_ 










* 












X 


' 


x 




CM 








a J 


*— N 


s 




Z 


x no m 


4 


a 


* 


O 










o 










CM 3 










CM l 


CM 


X 








*— • 


X 


CO 


X 




< 


1 


< 










¥ 






CM 






o 






X 


CO 






in 


o 


c 


< 




CO 


«*■ 


o 


1 

i 


o 


X 


3 


X 




t— 4 




tH 


u 


CO 




3 




no o 


X 


X ^ 


2 CD 


X 


Q0 


o 


4 


4 




1 


X 


X 


-k| 


tH 




CM 






X 


•—4 


X 




X 




rvi 








CO 


2 




•— * 


2 


; 4 


• 




< 


< 






H 




<=> 


•— * 


4 




4 


w | 


X 




a 


4 


tH 


CO 






*— * 


i 


o 


— » CO 


CO 




* o 


tH 


3 








*~s. ' 


►— » 


■k 




X 




x! 






CO 


•—4 


x. 




X 


\ 


i 


^4* 


CO 


o 


CO 




CO 


<c 


X 


X 


X 


X 


HH 


CM 


X 


c 


ro x o 


X 


O 


X 


CJ 


X 


CO 


CM X 


tH 




ro 


* 


• 


\ 




4 


4 






X 


1—4 


CD CO 


X 


w 


X 


CO X 


4 


tH 


4> CM 


X! 


4 


X 


¥ 


X 


z 


xi 


tH 


-j 


-~N rl 


x ro 


1 


I 




r 


l 




2 


4 


¥ 




— 


*. 


o 


CM 


a 


-J 


O 


— J 


— tH 




tH x 


CO 


a.1 


3 


;vj 




w* 


z -j 




X 


/*“V 


i 


1 


i— 


O 


ao 


CD 


X 


X 


Q 


X 




Q 


rvj 


Q 


rsi 


P: 1 


X 


X 


x; 


CO 


4 


CO, M 




N-3 


X 


- w 


X 




tH 


x— «» 


X 


CO 


3 







3 X 
X C i 

o 



o < 4 «< + 

• a: q a: a c\j o 
ro o I o 



Wa.2Q.UOO,^LffihW •> ^ x-i + -* _s 

3 — » x^ . Li_ cn z x w cm — 4 ^ x ^ <t o o 

COU.ll ll (/) O ^ O Ll X CD CZ CD 



o 




O X 


CM X 


1 II 


< 


II 


< 


tH 






• — 






CD 


O 3 


3 


X 






CD 


i 




X 


II 


II 




• 


<c 


< 






*— «. -i 


3 






X 




X 


II II 


II 


ll 


ll X 


4**kk 


#— «k 


< 






c 


X 


* 4 


l 


c 


3 


1—4 






UJ 


a 


X 


X 


r~ 




CO M" 


XI 


r— 


: '"*■ 








1-4 






3 


— « 


— • 


1 


ll II 


1 ll 


T— 


CO 


1 


< 






3 


X 


X 


X UJ 


< 


c 


< 




X 


* 


< 


tH 


II 


CM 11 


tH 


CM 


tH 


CM 3 






II 


j 


| 


<t o 


X 


kk^ 






N— ^ 




X 


X 


3 


II 


II 


X 




UJ Q 


X 


X 


’W' 








tH X 


X 


X 


X 


on 


3: 


3 X 


X 




o 




w 




ll 


1 

I 


X 


X 




3 


tH 


CM 


cc 




X < 


CM 


X 


<— « 


tH 


*— 4 


CM 


X 


X 


X 


X II 


X 3 rO 3 3 


X 


11 


**** 


1 11 


ll 


ii 




11 


3 


3 


CSC 




—4 


*— 4 


o 


X 


— LII 


3 


o 


. X 


3 


X 


3 


O CO 


CO X 


X 




o 


, 3 


O O 


4 

1 




X 








o o 


►—4 


•—4 


o 


UJ 


3 


3 


X 




Q X 


X 


X 


X 


XJ X 


XJ 


o o 


O 3 


3 3 


3 


o 


XJ 


o o 


3 


X 


*— t 


X 


p 

0 




00 


Q 


a x 


00 


X 


X 




CM 


rH 




o 
































! 





















in 



cm ro 



o o 



105 



40 I F ( ZL2 - E*l)7fl870Q6BO 
50 1 F ( ZL3 - PI ) 60,70.70 

60 „Z L 1 _=_Z L 2 

ZL2 = RADD*ZL1D + RADP *ZL1M ♦ RADS*ZL1S 
C0S3 ( 1 ) = C 0 S F ( P S I 2 ) 



(\1 H r( 

*— ♦ *— • *— < 

co in in 

cl Q. CL 
^ 

LlJ Ll_ U. 

2 co 2 

~ O — 

CO O OO 

II II II 

I £ 

<*->■ *-> rv 

H OJ CM 

O 

cnanonH 

2 W 2 ! 

*-• o — • o 

CO O 00 CD 



CL CL CL CL 

V* W 

L. L II L 
(^ 2(^2 
O — O — 
O CO O CO 
II II II II 







CO 


































* 


































X 












X 






















in 












2: 












•< 


1 








o 












►— » 












X 










o 












CO 












X 


































CO 








i 












s 












2 






















*-> 












-J 


1 
































a 


i 








CM 












T-l 












•i 










w 




1 




















X 








CO 












CD 












2 










2 












2 T 




*— * 








CD 




















1—4 












2 








OO 












OO 












-J 


1 








-w» 












* 




T-l 








a 










>-* 












X 












•< 






















CO 




CD — 








2 d 










1 












o 




OO CL 
























> 






o 




O \ 








> — 










• 






►— # 










O • 
















*H 






O 






1 




■H o 
























2 










X CO 








a 








* 
















2 . T-l 








•« 






_ 




X-S 






* 






CM 




— + 












oo 




CM 






CD 






w 




tn — 








2 






o 


• 


* 






1 






CD 












H 


i 




o 




+ 






LU 






2 




s ^ 








CO 




* 




*** 






O 






►— • 




<t 








•— ( 


i 




q 




s 






• 






OO 




^ Ll. 








o 




CM 




rl 






to 








! 


^ — 1 








1 


1 




w 




w 
















tH < 














CO 


\ 


CD 






% 










w | — 








1 






c n 




CO 






UJ 






1 




CD ^ 














q 




O 


1 


1 


CO 




/ 






2 Ll. 








j 






o 




o 


! 


1 


CL 






• 




— 2 








* 






♦ 




W 


1 




»—* 






T-l 


I 


tn <c 








Q 


[ 




H 










_J 






W* 




* i — 








CO 


\ 




tH 




* 






LU 






Ll 




X -Ii 








a 


| A 




w 


CM 




CM 


I CM 








I — 


! 


oo ~ 








Ql 




| 


no 


* 


CD 


*; 


* 


% 


X 




cr 




O T-l 




CM 






i o 


A 


oo 


♦ 


* 


♦ 


* 


X 


X 




o 




. c_> 




tH 




o 


tH 


•s 


o X 


CD 




*-s 


LU 


3 




co a 


1 GO 








CO 


Ll. 


O 


o 


CD 


1 


*— > 


<* 




CO 


t 




V 


^ 2 




r^«) 




CJ 


CM 


•L 


+ o 


<C 






*. 




• CM 




« 


CM — 




H 




0 -! 


« 


A 


*—n 


O 


+ 


tH 


1 Q 


• 


Q 


i A 


t-i 


a 


> 00 




•» 


• 


•> 


CM 


} 


CM 


V- * 




w 


O CM 


il 


00 




co 00 ^ 




ro 


o 


QD 


• 




H 


1 


W-' 


CD 


i X 


X 




• 


CD 


t— 


1 2 Ll. 




tH CO 


CO 


CD 




CD 




•f 


CO 


i CD 


»— i 


i 


T-l CO 




— 2 






tH O 


tH CM 


5 




CD 


o 




CL 


s 


O 




00 < 


a 






CL 


Ll. 






* 


* 


o 




LU 


X o 




w h- 


s 


O' 


1 






*— 4 


CO 


T-l 03 


* 


I 


o 


♦ 






CM ^ <t 


• 


-i 








.X 


* 


-s— ' 


w 


X 


/ 




< 

. 1 


o u. 


* 


Ll + 


o rsj 


ro 


< 


•k 


CL 



T-l I— 

^ DC 
CD O 
2 00 



tH 2 * 

>-» O 4 

cr i/o u. co fj 
►- q. ; 



tH tH CM CM CO II 

w w || 

no m no oo xi x 

C/) 2 (/) 2 (/) 2 

O — 1 O O 

O 00 O CO O 00 



a 

IS 



O wj 

oo s 

n ^ oo 

it 



^ ar x: it, 

^ G *—* i 

~ 00 oo s 



o 

O CD 

r 2 

CD >-? 



_J CD 
-J 2 ! 

c _i 

Q O Q 



2 co * oo tH go x co ro oo tH 

O ~ O -J tH I < O 2 < I 

Q U ^ U N ♦ J(l'< t LA 

<^ < v HfOX I X ^ 

II tH _J X * X O* tH <- 

ii w ii ii isi < ii co * ^ — ro 

x qq ; -J I < i- - J 

2 <t co <ro ii x fo i — x < in 

OXOL t | ^ ^ 2 3C 2 aj 
2 X O -J <t O' < — < 

-j <t * < j jl jcro 

Q CD X N — X CL 



CL 

O L ! U. Z 



i 

CM ro 



U. — « 

co 



Q 
2 
— LLI 

O' A 



Q 

2: 

•* 

cr 

CL 

ID 

O 

2: 

id 

Ll. 

•I 

00 

< 

CD 

* 

< 



o 

CO 

CL 



OO 



•JJ 

2 



2) 

O 

on 

CD 

to 

oo 



10 6 



CD O CD 
CD CD O 



CD CD 
O CD 
CD CD 



or 

o 

1X1 

z 



° 

III' 

3 
— I 
C 
> 

lu »— i— j or 
_J — — o 
3 Z Z 

or — — or 

3 3 ID 

co 3 i— 

z rt fr c /5 ?/) 

0 UU UJ Z UJ 

CO 3: CL. < I — 

non 

1 j: 

►-•ii it z 



o 

z 

< 

or 

o 

LU 



CO 

£ 

I _ 

Lu 

LU 

Q 



O 

j ►-• 

5 

CO LU 

z z 

O *-* *■ 

*— * 1 — 

CO 3 

— a 

> or 

— • CD 
Q 3 

J" 

O 7 ? 

! o 
or + 

LU — 

CD O 
21 Z 
3 3 
2 L H 
II II ll li 



O CD CD CD CD O 
CD O CD CD CD CD 
O O O CD O CD 



CM 

\ 

CD 

O 

z 

3 

Lu 

+ 

<£ 

0 
z 
3 
Lu 

* > 

<C Q 

1 Z 
QD * 
^ C\i 
ii n 



i 

u. c\j — x 



OU.O 
O — • 1 — 

3 Q <t 3 

Lu Z O tu 

V. * 3 + 

Lu O 
C ll * 3 

• I • 5 I — 

CD CO LU CO 

I|} w o □ II 



o o o o a u Q 

O 0 ,0 o o o o 
o oo o o o o 

O CD CD CD CD CD CD 



to 

* 



or 

3 

c_> 



ro 



> 

LU 

or 

CL 

I 

or 

3 

o / 
-/ 
lL 
:o 
CD 
<£ 






or 

3 

O 

+ 






UJ 

o 

* Lu 
• CO 

cd 

+ *< z: 
LU LU * or o 
3 z or 3 co 
z o or o ^ 

• CO LU II ’ ll 



ro 

CD \ 

tw or 
3 z 
o o cr 



LU 

co 

CL 



CO 

* 

XT 

CO 



LU 



or 

o 

co 



CM 

* 

* 



Q 



LU Z 

o co cr 



II II > 


o :> *> 


UJ ■> o 


ll CM + 


O »- 


II w t> 


Ul »— II 


3 O Z CL 3 


co or — 


iZ-LU 


12-2 


3 < 


3 z or lu z co 


1 - o z 


2 : t* } — 


z or a 


i 3 o or 


O O *— LU O ll 


1 — O 3 Lu or 


O O Z UJ Z 3 


O — l LU 


K 1 UZ 


; Ll. Z CL 


CO Z CO 


Q Q X 


CO o 


O — CL 


CO o c 

j 


a: LU Lu 


o lu or 






o 


i 


CD 


' CD 


<=> 




! 






T-i 


1 

1 


CM 


fO 








* o o o 


O 






* 












i 

i 






M 

* 








j 



Q Q 

z 

LU LU 



107 



>- 

UJ 

'X 



x 

>- UJ 
UJ h- 
X 2 T 
UJ o 
H- X 

o o 



< 

x 

o 

-J 

< 



o 

< 

a 



x 
>- o 

LU 51 

x 

UJ O 
h- h- 

2: 

O h- 
2; x 
o 

O Cl 



UJ 

x 

UJ 



>- h- 

2: 
o 



x 



X 



m o 

D U I- 
^ O' 

~ x x 

< O CO 

X u_ o 



> 

UJ >- > > 

DC UJ UJ UJ 

UJ X X X 

UJ 2T h~ UJ UJ LU 

X I — I — H- 

O X X X 

5 ! O O O 

51 SI SI 

o 

o a h- o o o 

X h- h- H- 



2: 

o 






< > to < 2: 

< J S in CD CL 

x x :s x x x 



00000000000 

cncncncncncncncncncnco 

cmcmcmcmcmcmcmajcmcmcm 

mLnmmmintnmmmin 



CSJ f\J f\J CSJ f\J f\J f\J CSI fNJ CM OJ 



OvOnOOvOsO^vO^^vO 

cncncncncncncncncncncn 

vDvOOvDOOOOOOO 

cncncncnmcncncncncncn 



l 

0 
o 

1 

m 



O CO o 

CM <J* r-l 



CM 



coomn^mcnr^chO'' 
co ^ o h h cm m 

r-tf^mr-lf^-r-imvOChCO 

omr-ovOCMOf^*r^iA 



D O H H 
H O f\J <f 



m 



incvjstvo. c\jO'fvjoo'>fiA 
CM<rCMCMCMCOH<riA CM 

vOo»— icncsjvfcoocoo^H 

vOr-tC\J^iT\vf 4 <f(f| (M 



3 

P 4 

fl 



© 

rH 

0 , 

§ 

W 









CM CO 


ro 


As 


tH 


1 i-f 


ro 


^ ro 


ao 


LA 










As 


LA 


ro 


CM 


YT 


co 


H 




O 








y-» As 


LA 


O 




On 


O 


V 








LA O 


0 


ro 


yH 


1 CM 


r-H 


y-t A* 


Y“i 


vO 


1 


3 






0 ro 


'r-f 


0 


NT 


M- 


As 


LA 


CD 


CO 


rH 


1 • 


K 






CM ro 


O 


vO 


t-H 


LA 


00 


00 


NO 


YT 


nO 


< 


c 




ro to 


ro 


ro 


ro 


ro 


ro 


ro 


ro 


ro 


ro 


a 


















l 






i 


t 


0 




O 


r*v 


On 


On 


GO 


CM 


ro 


r'v 


LA 


rH 


1 


X X 


co 0 


0 




O 


O 


CO 


ro 


H 


rH 


LA 


a 


1 


























_ 


J 


k- 


As 0 


tH LA 


ro 


co 


NO 


>o 


rH 


NO 


As 








3 


0 CM 




LA 


NO 


CM 


O 


ON 


CO 


O 


CM 






c 


ro 0 in 


yT 


On 


Y-» 


O 




NT 


As 


LA 






X 


As 00 


O 


On 


NO 


On 


O 


0 


On 


00 


On 








rH h 

I 




t~I 


iH 


<r-l 


CM 


CM 

i 


tH 


rH 


yH 


1 








1 

a* in 


CD 


M" 


As 


As 


LA 


1 

iH 


LA 


CD 


O 










a* co 


LA 




ro 


CM 


in 


YT 


CM 


ro 


M" 








X N CO 


O 


T“H 


ro 


NO 


00 


CD 


GO 


ro 


ro 








LTV LA LA 


H 


CD 


t-t 


CO 


ro 


NO 


ro 


ao 










O O 


CO 


CD 


O 


NO 


CM 


LA 


On 


CM 


CO 


0 




O 


ro 


yt 


00 


yJ- 


yH 


ro 

1 


LA 


ro 




CO 




i 












l 










UJ 




i 






















C 


\ 




LA O 




f"s 


LA 


00 


NO 


LA 


GO 


LA 


As 




0 X 


00 ro 


t-i 


M- 


LA 


As 


ro 


NO 




co 




LU 


























O ►— 


nO LA 


CM 


NO 


in 


CM 


IA 


NO 


CM 


LA 


vO 




\ 


3 


O O .00 


NO 


yT 


On 


On 


CM 


YT 


CM 


On 






c 


CM IA 


CD 


tH 


vO 


5 T 


NO 


co 


H 


CO 


O 








X 


m- ro 


CM 


CM 


M" 


CM 






CM 


CM 


CM 



/ 



I 

>-, 

LU 

or- 

>- LU >- 
LU J— UJ 

or 3 or 

UJ O LU 



OOO 

xi — 2: 



>- 

LU 

QL 

LU 

}— 

7 

O 

X 



>- 

UJ 

cr 

UJ 



> 

LU 

X 

LU 



< 

X 



OOO 
h- c h- 

o 

<c 

cr x { 

□ ^ x or 

J X < O CD 
Cl — X Ll C5 



x o 

O H~ 
X 



CO 

UJ 

cr 



o o 



< M 
C _J 



i 



LU >- > >- 
rr lu uj lu 
uj cr cc cr 

I — LU LLI LU 

X * { — ! ^ 

O X X X 

x o o o 
XXX 

o 

1-000 



> cn < x 

X CO CD CL 
T X X X 



3 

CL 

•P 

3 

0 

© 

iH 

Ot 

1 

CO 



108 



APPENDIX II 



PROPORTION OF DAYLIGHT AND DARKNESS ON PATHS OF PROPAGATION 

The difference in the height of the D-layer of the ionosphere from 
day to night is the largest single cause of variations in the phase vel- 
ocities of VLF waves and causes variations in the attenuation rates of 
a magnitude equal to those caused by other factors. Because of this the 
percentage of daylight and darkness on each path of interest must be 
calculated as a function of time. 

This problem has been solved by Brady and Crombie £223* However, 
since the following solution was developed independently, it will be out- 
lined here. 

The following is based on the premise that the percentage of day- 
light on a path can be found if four points are known: (1) transmitter 
location, (2) receiver location, (3) and (4) intersection of the sunrise 
and sunset lines with the great circle path between (1) and (2). (A 
spherical earth of unspecified radius is assumed.) 

The great circle path can be described by the longitude of one of 
the points where it crosses the equator (A^) and its declination at 
the equator . The sunrise-sunset line forms another great circle 

that can be described in the same way. The intersection of these two 
great circles on the sunrise side of the earth will determine the coor- 
dinates of point (3) mentioned previously. The geometry of the situation 
is shown in Figure II-l. 

Here: 

A 2 = 15t ~ 90 degrees 
OC 2 = <x s + 90 degrees 
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where 



t 






= GMT 

= declination of the sun 




Using spherical trigonometry the coordinates of P (^3? X 3) can 
easily be found: 



■1 r . cosigCtxj - eK 2 ) 

a = tan tan%(A 2 “ A p 

L cos%( + X 2 ) 



'[ 



+ tan 1 tan%( \ 



sin^CfrCx 
1 sin%(K j 



-*2> I 

+ * 2 > J 



tan c = tan a cos(180° - ~ - tan a cos (X 2 

^ 3 = A 2 + c 

sin j 2 3 = sin a sin 0< 2 

The point at which the sunset line crosses the great circle path 
is then: 

A 4 = X 3 1 180° 

$4 = -&Z 



At a height of 70 km the point P will be further along the great 
circle path in a direction away from the sun. Figure II-2 shows the 
arc distance by which the sunrise and sunset lines are shifted. 

Here: 

h = ionospheric height 
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S = screening height 



R = earth radius 





sun s rays 



Figure II-2. Sunrise/sunset at ionospheric height 



The screening height, mentioned by Brady and Crombie £ 22 ]}, is the 
height below which the components of the sun's rays that most affect 
ionization cannot penetrate. This, in effect, changes the radius of 
the earth by an amount S. From Figure II-2 then: 



The point at which this new sunrise line crosses the great circle 
path (P^) is shown in Figure II-3. The coordinates of this point can 
be found by spherical trigonometry: 



tan%(a + d) = tarfg(X 2 “* Xp 



cos%( i - *2^ 



cos^C^i + 2 ^ 



tan%(a - d) = tan %(^2 " Xl^ 



sin^(OC i - 



sin%( K i + b< 2 ) 



coth f = tan %(^< 1 - &< 2 ) 



sin%(a + d) 



sin%(a - d) 
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Figure II-3. Revised position of sunrise point 



From Figure II -4: 

cos = sin Xj cos( - \ j) 

^ ^ cot Ad cos fiz - s in #3 cos /° 

sin^O 

A3 = A3 + 0 



sin 



/3 3 = sin 3 cos Ad + cos jS ^ sinAd cos p 



P0L6 




The coordinates of the point of sunset A4) can b e found in 

a similar manner from Figures II-5 and II-6: 

X4 = As i 18 o° 
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Ad 
cot 9 ' 



= -Ad 



cot Ad' cos ^ + sin cos 
sin^A 



A 4 = A 4 + 0/ 



sinfi^ = sin ft ^ cosAd - cosy ^4 sinAd cos^O 





Figure II-6. Sunset position coordinates 
It should be noted here that if it is assumed that (2 ^ and Ad ' 
are negative quantities, then 9' will also be negative and the resulting 
longitude correction will be in the opposite direction from the sunrise 
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longitude correction. 

Once the coordinates of the four required points are determined 
the amount of sunlight can easily be found with the aid of Figure II-7. 
This shows the four possible situations that can exist on the great 
circle path of propagation. Points (1) and (2) are the transmitter and 
receiver positions, (3) is the sunrise point, and (4) is the sunset 
point. The amount of daylight on the short path for case A is the dis- 
tance from (3) to (2); the amount of daylight on the long path for case 
A is the distance from (4) to (2), and so on. 



om. 

case D 

Figure II-7. Possible positions of sunrise and sunset 
on great circle path 

The actual program is shown on pages 116 through 120. The geodesic 
program was converted to a subroutine for convenience here. In the pro- 
gram the declination of the sun as a function of the time of the year 
is approximated by two sine waves of different periods; one being used 
from the 1st to the 80th and the 267th to the 365th days, and the other 
being used from the 81st to the 266th day. Data from the Nautical Al- 
manac C173 was used as a guide in setting up this approximation. 

The program also includes a correction factor for the Greenwich 
Hour Angle, which was also found necessary from data in the Nautical 
Almanac. This correction factor has a maximum value of 4.1 degrees 
(this changes ^2 i- n Figure II-l) and is again approximated by sine 





case A case B case C 
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waves of various periods and amplitudes (see statements 311, 511, 711, 
and 812 in the program). 

Geocentric latitudes, obtained from the geodesic program, are used 
so that no error is involved due to the ellipticity of the earth. The 
radius of the earth is not even used in this program, since the distances 
between points are expressed as arc distances and the amount of daylight 
is expressed as a percentage of the long or short path. 

The only inputs required are the time of day and the day of the 
year. The latitude and longitude of the two points of interest and 
the description of the great circle path through these points are ob- 
tained from the geodesic subroutine. 

In the paper by Brady and CrombieC22^ three basic simplifications 
were made which could result in some error. These were the assumptions 
that: 

1) the sun is a point source, 

2) atmospheric refraction can be neglected, 

3) the earth is spherical. 

The third assumption does not introduce any error in the program devel- 
oped here because geocentric latitudes were used. The other two could, 
according to Brady and Crombie, result in errors of about 16 minutes 
and 34 minutes of arc respectively in the quantity Ad (Figure II-3). 
These errors, although small, could cause considerable error on high 
latitude great circle paths. 

This program is not valid for the case in which the highest latitude 
on the path is greater than the quantity (90° - bC g ). This causes a por- 
tion of the path to be continually in daylight or darkness^, (e.g., polar 
paths during the winter or summer) and the method used here breaks down. 
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JOB0564F, MCKAY, J.D.EAY4 
PROGRAM SHADOW 

this program computes the percentage of darkness and daylight on 

The LONG AND SHORT PATHS (GREAT CIRCLES) BETWEEN TWO POINTS ON 
THE SURFACE OF THE EARTH FOR ANY TIME OF DAY ( UT ) AND ANY DAY OF 
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1 + PI/2. 

GO TO 100 

90 ALFA ( 2 ) =SINF( ( 360 . /357 . 2 > * < D ( I ) -8 7. 9 ) *RAD ) * ( 23 . + ( 2 6 . 7/60 . ) ) *RAD 
1 + PI/2. 
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FOLLOWING IS A POSITION CORRECTION TO OBTAIN THE SR AND SS 
POSITIONS AT A HEIGHT OF 70 KM ABOVE EARTH. 
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APPENDIX III 



LOOP ANTENNA CALCULATIONS 



Consider Figure III-l below. This is a schematic representation 
of a loop antenna showing a voltage generator in each vertical side of 
the loop to indicate the voltage induced by the passing wave. The z- 
axis is taken as the vertical axis and the xy plane the horizontal plane 
with respect to the earth. The voltage E^ appearing across the tuning 




Figure III-l. Loop Antenna 

capacitor C is the signal produced by the antenna and is found by phasor 
addition to be E^ = E 2 - E^. The symbol ~ here denotes a phasor quan- 
tity, while a letter without this indicates magnitude only. 

In order to find the effect on the antenna output of various di- 
rections of arrival of the signal, assume the antenna to be rotatable 
about the z-axis. The angle (X is defined as the angle between the 
normal to the loop and the direction of arrival of the signal, as shown 




Dl£E.CT\OM OF 



ARRIVAL 



Figure III -2. Loop Geometry 
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The angle varies from 0° to 360° and taken to increase positively as the 
loop rotates clockwise. The bar shown on the loop in Figure III-2 is 
for reference and corresponds to side 2 of the loop as seen in Figure 



III-l. 

The phasor diagram of the loop voltages is given in Figure III-3. 
The angle (/? represents the phase of the antenna output signal, where 
4? = 0 when the loop is oriented broadside to the direction of arrival 
of the wavefront. 




Figure III -3. Phasor Diagram 

The angle 0 in Figure III-3 represents the difference in phase be- 
tween the voltage induced in side 2 and that in side 1 of the loop. It 
is determined by the distance d between the two vertical sides, and is 
given by 

0 = - - 71 — ■ sinOt radians, where 

A 

A = free-space wavelength of the received signal 
(same units as d). 

Since = E 2 > the phasor diagram forms an isosceles triangle and there- 
fore = ^9. The phase angle is then 

C P - ■ j. sinOC radians. 
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The resultant is also found readily from the phasor diagram. 
Using the Law of Cosines, one has 
2 9 2 

= + E2 * 2E^E2 cos ® 

= 2Ei - 2E^ cos 9 
= 2E^ (1 - cos 0) 

E r = E l H (1 - cos 0) 

= Ei V2 (VT sin 



= 2E^ sin (p . 

It is convenient to normalize by setting E^ = 1. One can then express 
the magnitude of the relative output E QUt in db below E^ ^ : 



E out 



20 



l°gl0 



e R 

E Rmax 



= 20 log. n ( | S f n % ), where If 

10 2 sin ip Q TO 






= 20 logjg 



. sin (P 
sin (f c 



■) • 
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APPENDIX IV 



CARDIOID ANTENNA CALCULATIONS 

As in Appendix III, assume the direction of travel of the incoming 
signal fixed and the loop antenna rotatable, with the angle OC defined 
as before; see Figure IV- 1. In addition, the sense antenna is placed 
at some arbitrary position, forming the angle tiC as shown. The distance 
d s of the sense antenna from the loop midpoint is non-crit ical , so long 
as it is short with respect to a wavelength. 




Figure IV- 1. Card io id Geometry. 

Because of the 180° phase reversal of the loop antenna when 
(X = 0° and 180°, it is convenient to analyze the cardioid system in 
two parts, with 0° < OL < 180° and with 180°< 0C < 360° (or equivalently, 
-180°< 0 i < 0°) . As in Appendix III, the loop phasor resultant E R is 
given by E^ - E^ and phase ip is 0/2. This is illustrated in the phasor 
diagram of Figure IV-2. In this diagram, both 0 and vf are taken to be 
negative angles. 

For proper adjustment of the cardioid pattern, the magnitude and 
phase of the signal from the sense antenna must be made equal, respec- 
tively, to the magnitude and phase of the signal from the loop antenna 
when the loop is positioned for maximum response ( Oi = 90°). This 
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Figure IV-2. Loop Phasor Diagram 

condition is shown in Figure IV-3, where E g is the magnitude (after the 

necessary amplification and attenuation) of the sense antenna signal. 

The angle (p o is the largest (negative) value of tp . Once set, the 

/\> 

magnitude and phase of E g , of course, remain constant as the loop is 
rotated . 

fr 

X 

E* 

Figure IV-3. Cardioid Adjustment 

The phasor diagram for some value of OC between 0° and 180°, other 
than <X = 90°, is shown below in Figure IV-4. The overall signal out- 
put of the cardioid unit is given by E rt /f° (note that T is negative and 
its magnitude less than that of CP Q ). This is obtained by the phasor 
relationship E^ = + E g , accomplished in the cardioid unit which mixes 

the loop and sense antenna signals. 
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X 




Figure IV-4. Card io id Phasor Diagram 
As in Appendix III, the angles 4?, <P o , 0, and 0 Q may now be ex- 
pressed in terms of the angle of loop orientation OC , as follows: 



9 



0 



o 



( D 
Y o 



__ -d sin CX 

A 



-d sin CL 

A 



max 




= 



max 



-d 

A 



x 360° 
x 180° 
x 360° 

- x 1800 



(degrees) 



The objective is to find 'Y and E q in terms of CX . At this point, 
it is convenient to derive an expression for Er as a function of E g . 
This is done as follows, referring to Figure IV-2 and using the Law of 
Cos ines : 

E| = E^ + E^ - 2E 1 E 2 cos 0 

= 2e| (1 - cos 9), since Ej = E2. 

E^ = E^ V7 s in 0 



0 = sin 0, 




© vr. 



since 0 is small. 
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ERma * _ Q ~ ES 

Ej o ' z ~ E 1 ' 

e r _ _e_ 

e s 0 O 



Figure IV-5 is an expanded view of the phasor diagram of Figure 
IV-4, drawn on an exaggerated scale to make angular relationships 
clearer* In this diagram and the following equations, the fact that 




and are negative is taken into account. In the triangle of Fig- 
ure IV-5, the two sides E R and E g and the angle between them are known, 
so that ^ may be found by geometry. One can then determine *T = ^ ^ , 
which is the phase angle of E q . Solving the triangle, 



^ = 90° 



180° + tp n - <£ 

2 



+ tan 



Es ~ e R 
E s + E R 



cot( 1800 'fn-V > 

= 90° - 90° + — ~ - + tan -1 



cot (90° - 



4> - (p 



127 




1 - X 



1 + X 




where x = -q— 
y o 



Since and ^ can be found for any given value of CX , T is uniquely 
determined. It is a negative angle and cp Q <T < 0° for 0° < & < 180°. 

Determining the relative amplitude response for these values of 0C 
is very straightforward. Refer to Figure IV-5 and use the Law of Cosines 
as follows: 



One may now normalize by setting E s = 1 and taking 201og 1 QE Q to obtain 
relative output in db, with zero db occurring at 0L =0° and 180°: 



The procedures for obtaining cardioid phase and amplitude for loop 
positions 180° < 0C < 360° are quite similar to the methods given above. 
Figures IV-6 and IV-7 below are self-explanatory and are analogous to 
Figures IV-2 and IV-4, respectively, discussed above. 

Angles are similar to those determined before, the principal dif- 
ference being that 180° is now added to ^ due to the difference in loop 
position. Note that all of the following angles are now positive, while 
y remains negative: 



Eg = e| + E 2 - 2E r E s cos [l80° + ( (f Q - (f )] 
= E| + E 2 + 2E R E S cos (Cf? 0 - 4? ) 




cos ( if Q - ^ ) 






db. 
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Figure IV-6. Loop Phasors (180° < CL < 360°) 




Figure IV-7. Cardioid Phasors (180° < (X < 360°) 



9 = - sin 0. x 360° 

A 



= 180° - — sin a x 180° 

A 



_d_ 

A 

d 



(£ 0 = x 180° 



0 O = x 360° . 



As before, ^ is found by trigonometry as follows: 



_ n tp n + if - 180°) _i 

X = 90° - — S-2 1 + tan L 

’ 2 



cot Vo * V . 180° 



E S - E R 

E s + e r 



= 180° - — - + tan -1 



x-1 4>o + 

— tan ^ 
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If I = 180° - [( q? - 180°) + 

T = - |l80° - [( - 180°) + 

= -180° + If - 180° + 

T = ip + \ - 360°. 

Relative amplitude is found as before using the Law of Cosines, and 
referring to Figure IV-7 : 

E l= E l +E l~ 2 E R E S c ° s [ 180 ° - < + t * >] 

^°- ■ = V 1 + X 2 + 2x cos (tf Q + tf ) . 
s 

Normalizing and expressing the result in db, one obtains: 

E q = 20 log^o 'si 1 + x 2 + 2x cos(tf Q + If ) db. 
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